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Abstract. We derive a systematic high-frequency expansion for the effective 
Hamiltonian and the micromotion operator of periodically driven quantum 
systems. Our approach is based on the block diagonalization of the quasienergy 
operator in the extended Floquet Hilbert space by means of degenerate 
perturbation theory. The final results are equivalent to those obtained within 
a different approach [Phys. Rev. A 68, 013820 (2003), Phys. Rev. X 4, 031027 
(2014)] and can also be related to the Floquet-Magnus expansion [J. Phys. A 
34, 3379 (2000)]. We discuss that the dependence on the driving phase, which 
plagues the latter, can lead to artifactual symmetry breaking. The high-frequency 
approach is illustrated using the example of a periodically driven Hubbard model. 
Moreover, we discuss the nature of the approximation and its limitations for 
systems of many interacting particles. 
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In the last years the concept of Floquet engineering has gained more and more interest. 
This form of quantum engineering is based on the fact that the time evolution of a 
periodically driven quantum system is, apart from a micromotion described by a time- 
periodic unitary operator, governed by a time-independent effective Hamiltonian l][2 


The aim is to engineer the properties of the effective Hamiltonian by designing 
a suitable time-periodic driving protocol. This concept has been employed very 
successfully in various experiments with ultracold atoms in driven optical lattices. This 
includes dynamic localization [S IO , “photon”-assisted tunneling 11-18 , the control 


of the bosonic superfluid-to-Mott-insulator transition 
Bloch bands 


19 20 


21-24 , the dynamic creation of kinetic frustration 25 26 


resonant coupling of 
as well as 


the realization of artificial magnetic fields and topological band structures 26 - ^ (see 
also Ref. for the creation of a topological band structure in an array of optical 
wave guides). In a quantum gas without a lattice, periodic driving has recently also 
been employed to tune or induce spin orbit coupling. 

A prerequisite for Floquet engineering is a theoretical method to compute the 
effective Hamiltonian (as well as the micromotion operator), at least within a suitable 
approximation. In the high-frequency limit a rotating-wave-type approximation can 
be employed for this purpose. This approximation coincides with the leading order 
of a systematic high-frequency expansion that provides also higher-order corrections 
to the effective Hamiltonian and the micromotion operator 42 -45 . In this paper we 


show that this high-frequency expansion can be obtained by employing degenerate 
perturbation theory in the extended Floquet Hilbert space. Our approach provides an 
intuitive picture of the nature of the approximation and the conditions under which 
it can be expected to provide a suitable description of a driven quantum system. We 
point out that the time scale on which the approximation is valid can be increased 
by increasing the order of the approximation for the effective Hamiltonian, while 
keeping a lower-order approximation for the time-periodic micromotion operator. We 
also address the relation between the high-frequency expansion derived here and the 
Floquet-Magnus expansion (see also Refs. 47-^). The origin of a spurious 


dependence of the quasienergy spectrum in Floquet-Magnus approximation on the 
driving phase is discussed (see also references 42 ^ ^ ). Using the example of a 
circularly driven tight-binding lattice, this artifact is, moreover, shown to produce 
a non-physical breaking of the rotational symmetry in the approximate quasienergy 
band structure. Finally, we discuss the validity of the high-frequency approximation 
for systems of many interacting particles. 

This paper is organized as follows. Section gives a brief introduction to the 
theory of periodically driven quantum systems (Floquet theory) and serves to define 
our notation. In Sectionwe formulate the problem that is then attacked in Section]^ 
by means of the degenerate perturbation theory developed in Appendix C[ The relation 
to the Floquet-Magnus expansion is discussed in Section [5 and Section illustrates 
the approximation scheme using the example of a circularly driven hexagonal lattice 
Finally Section discusses effects of interactions within and beyond the 


27 38 39 


high-frequency approximation, before we close with a brief summary in Section 
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2. Quantum Floquet theory and notation 

2.1. Floquet states 

A quantum system described by a time-periodic Hamiltonian 

H{t) = Hit-^T) (1) 

possesses generalized stationary states |V^n(t)) called Floquet states [^. These states 
are solutions to the time-dependent Schrodinger equation 

indtlipit)) = Hit)\'ilj{t)) (2) 

of the form 

\Mt)) = |w„(t))e—(3) 
with real quasienergy and time-periodic Floquet mode 

\Un(t)) = \Un{t -\-T)). (4) 


Here dj denotes the derivative with respect to the time t. The existence of Floquet 
states in time-periodically driven systems follows from Floquet’s theorem in a similar 
way as the existence of Bloch states in spatially periodic systems. For completeness, 
we give a simple proof for the existence of Floquet states in [Appendix A| 

The Floquet states are eigenstates of the time-evolution operator over one driving 
period, 

U{to-\-T,to)\ipnito)) = (5) 

Here U{t 2 ,ti) denotes the time evolution operator from time ti to time ^ 2 - The 
eigenvalue jQgg depend on the time tp from which the evolution over 

one driving period starts. Therefore, one can obtain the quasienergy spectrum by 
computing and diagonalizing Uifa -\- T,to) for an arbitrary tg. The time-dependent 
Floquet states lifnit)) can subsequently be computed by applying the time-evolution 
operator, \tfn{t)} = U{t,to)\'tlJn{to))- 

The Floquet states can be chosen to form a complete orthonormal basis at any 
fixed time t. As a consequence, the time evolution operator can be written like 

n 

Moreover, one can express the time evolution of a state |'0(t)) like 

n 

with time-independent coefficients c„ = (un(to)|'0(io))- That is, if the system is 
prepared in a single Floquet state, |c„| = 5n,noi time evolution will be periodic 
and (apart from the irrelevant overall phase factor described by the Floquet 

mode (t)). If the system is prepared in a coherent superposition of several Floquet 
states, the time evolution will not be periodic anymore and be determined by two 
contributions. The first contribution stems from the periodic time dependence of the 
Floquet modes |u„(t)) and is called micromotion. The second contribution, which 
leads to deviations from a periodic evolution, originates from the relative dephasing 
of the factors Thus, beyond the periodic micromotion, the time evolution 

of a Floquet system is governed by the quasienergies e„ of the Floquet states in a 
similar way as the time evolution of an autonomous system (with time-independent 
Hamiltonian) is governed by the energies of the stationary states. 
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In order to study the dynamics over time spans that are long compared to a single 
driving period, one can ignore the micromotion by studying the time evolution in 
a stroboscopic fashion in steps of the driving period T. Such a stroboscopic time 
evolution is described by the time-independent Floquet Hamiltonian Htg- defined 

such that it generates the time evolution over one period, 

e^p(^-'^TH[)j=U{t^ + TM)- ( 8 ) 

and can be expressed like 

Hto=^^n\Uu{to)){UuitQ)\. (9) 

n 

The parametric dependence on the initial time to is periodic, = Hf^, and related 

to the micromotion. It indicates when during the driving period the dynamics sets 
in or is looked at and should not be confused with a time dependence of the Floquet 
Hamiltonian. From a Floquet Hamiltonian Hf^ obtained for the initial time to one can 
construct a Floquet Hamiltonian for a different initial time tg by applying a unitary 
transformation, H^ = {to,t'o)H[g(j{to,t'o). 

It is convenient to introduce a unitary operator that describes the periodic time 
dependence of the Floquet modes, i.e. the micromotion. Such a two-point micromotion 
operator can be defined by 

UF{t2,ti) = '^\Un{t2)){Un{tl)\ ( 10 ) 


SO that, by construction, it evolves the Floquet modes in time, 

\Un{t2)) =UF{t2,ti)\Un{tl)). ( 11 ) 

It is periodic in both arguments, C/f (^2 +T,ti) = UF{t 2 ,ti -\-T) = UF{t 2 ,ti). 

If the Floquet states and their quasienergies are known, e.g. from computing and 
diagonalizing the time evolution operator over one period, one can immediately write 
down the Floquet Hamiltonian and the micromotion operator by making use of Eqs. 0 
and (10). However, both the Floquet Hamiltonian Hf^ and the micromotion operator 
UF{t,t') might also be computed directly, without computing the Floquet states and 
the quasienergies before. This will be the aim of the approximation scheme described 
in the main part of this paper. From the Floquet Hamiltonian and the micromotion 
operator one can then immediately write down the time evolution operator like 

U{t2,ti) = H) = UF{t2, (12) 


Moreover, the Floquet modes |u„(to)) and their quasienergies can, in a subsequent 
step, be obtained from the diagonalization of Hf^^, 

HfjUnito)) = en\Un{to)). (13) 

The periodic time-dependence of the Floquet modes can subsequently be computed 
by employing the micromotion operator, \un(t)) = UF{t,to)\un{to)). 
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2.3. Quasienergy eigenvalue problem and extended Floquet Hilbert space 

The phase factors and the Floquet states |'(/'ra(t)), solving the eigenvalue 

problem of the time-evolution operator over one period, are uniquely defined (apart 
from the freedom to multiply each Floquet state by a time independent phase 
factor). In turn, the quasienergies e„, and with them also the Floquet modes 
\unit)) = and the Floquet Hamiltonian are not defined uniquely. 

Namely, adding an integer multiple of hw to the quasienergy £„ does not alter the phase 
factor Fixing each quasienergy £„ within this freedom fixes also the Floquet 

modes and the Floquet Hamiltonian. For example, one can choose all quasienergies 
to lie within the same interval of width hui, often called a Brillouin zone. This term 
reflects a loose analogy to the theory of spatially periodic Hamiltonians, where the 
quasimomentum can be chosen to lie within a single reciprocal lattice cell such as the 
first Brillouin zone. 

Starting from the known solution given by |u„(t)) and £„, one can label all possible 
choices for the quasienergy by introducing the integer index m, 

enm = £n + mhu!. (14) 


The corresponding Floquet mode reads 

\Unmit)} = 

such that 

\Mt)) = |«n(t))e-*^"*/'^ = |u„„^(t))e— 


(15) 


(16) 


When entering the right-hand side of Eq. (16) into the time-dependent 
Schrodinger equation ([^, we arrive at 

[H (t^ l^d.t\\Ujim(t)} = ^nm l'^nm(t)) ■ (1^) 


This equation constitutes an eigenvalue problem in an extended Hilbert space H = 
T~L ® C-T 01 This space is given by the product space of the state space H. of 
a quantum system and the space of square-integrable T-periodically time-dependent 
functions Time is treated as a coordinate under periodic boundary conditions. In 
the extended Floquet Hilbert space J-, the scalar product combines the scalar product 
of H with time averaging and is defined by 

{{u\v)) = ^ dt{u{t)\v{t)). (18) 

We will use a double ket notation |m)) for elements of H; the corresponding state at 
time t inH will be denoted by \u{t)). Vice versa, a state |u(t)) = \v{t-\-T)), including 
its full periodic time dependence, is denoted by |u)) when considered as element of H. 
In the following we will stick to this convention and conveniently switch between both 
representations. Likewise, an operator acting in J- will be indicated by an overbar to 
distinguish it from operators acting in "H, which are marked by a hat. For example, 
Q denotes the J^-space operator that in Fi is represented by 

Q{t) = H{t) - ihdt. (19) 

The operator Q is called quasienergy operator. It is hermitian in H and, as can be 
inferred from Eq. 0, its eigenstates and eigenvalues are the Floquet modes and their 
quasienergies, 
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Figure 1. Block structure of the quasienergy operator Q with respect to the 
“photon” index m. Each block corresponds to an operator Qm'm — ^m'—m + 
S^i^mhu) acting in the full state space H. The diagonal blocks Ho + mhu) can 
be interpreted to act in the subspace of relative “photon” number m and the 
off-diagonal blocks which obey H^i_^ = describe (m' — m)- 

“photon” processes. 


The complete set of solutions of the quasienergy eigenvalue problem (20) contains a 
lot of redundant information. In the extended space \unrrif) and |Mnm')) constitute 
independent orthogonal solutions if m' ^ m. These solutions are, however, related 
to each other by Eqs. (14) and (15), and give rise to the same Floquet state |'0„(t)). 
All Floquet states |V'ra(0) of the system can, thus, be constructed, e.g., from those 
Floquet modes whose quasienergies lie in a single Brillouin zone of the /ioj-periodic 
quasienergy spectrum. 

The quasienergy eigenvalue problem (20) 
computing the Floquet states or the Floquet 


provides a second approach for 
Hamiltonian, alternative to the 


computation and diagonalization of the time evolution operator over one driving 
period. It provides the Floquet modes not only at a time t^, but including their full 
periodic time dependence. Despite the drastically increased Hilbert space, treating the 
quasienergy eigenvalue problem (20) has also advantages. In order to diagonalize the 
hermitian quasienergy operator Q, one can employ methods, concepts, and intuition 
from the physics of systems with time-independent Hamiltonians. When describing 
parameter variations, such as a smooth switching on of the driving amplitude, one can 
even derive a Schrodinger-type evolution equation acting in Floquet space and apply 


the adiabatic principle 


X A discussion of smooth parameter variations in a driven many-body lattice system can be found 
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A complete set of orthonormal basis states \am'f) of H can be constructed by 
combining a complete set of orthonormal basis states |a) of "H with the complete set 
of time-periodic functions labeled by the integer m, 

\am{t)) = (21) 

From this restricted class of basis states all possible sets of basis states can now be 
constructed by applying unitary operators C/, \am'f)ij = U\amf). With respect to the 
basis jam)) the quasienergy operator possesses the matrix elements 

{{a'm'\Q\am)) = 1^ [ dt 

T Jo 

= (o; \H'rn' “t“ dfri'mda.'^ (22) 


where 


Hm — rp 


dte 


^H{t) = H\ 


(23) 


is the Fourier transform of the Hamiltonian H{t), such that H{t) = J2m=-oo 
With respect to the Fourier indices m the quasienergy operator possesses the 
transparent block structure depicted in Fig. Each block represents an operator 
Qm'm = Hm'-m + Sm'mmhw acting in T-L. 

The structure of the quasienergy operator Q resembles that of the Hamiltonian 
describing a quantum system with Hilbert space H coupled to a photon-like mode in 
the classical limit of large photon numbers, where the spectrum becomes periodic in 
energy. In this picture m plays the role of a relative photon number. The quasienergy 


eigenvalue problem (201 is, thus, closely related to the dressed-atom picture 52 53 


for a quantum system driven by coherent radiation 54 . Based on this analogy, one 
often uses the jargon to call m the “photon” number. Moreover, the matrix elements 
of Hm are said to describe m-“photon” processes. This terminology suggests a very 
intuitive picture for the physics of time-periodically driven quantum systems and is 
also employed when the system is actually not driven by a photon mode. 

In order to diagonalize or block diagonalize the quasienergy operator, it is natural 
and sufficient to consider unitary operators U that are translationally invariant with 
respect to the photon index m, {{a'm'\U\am)) = {a'\Um'-m\o)■ They correspond 
to time-periodic unitary operators U{t) = X]m=-oo acting in H (see also 


Appendix B). From Eq. (19) we can infer that a unitary transformation with such an 


operator [/, 


Q 

1^)) 


Q' = Wqu 
\u')) = W\u)), 


(24) 


is equivalent to a gauge transformation 

H{t) H\t) = U\t)H{t)U{t) - ihU'<{t)dttJit) 

= ( 25 ) 

with a time-periodic unitary operator U{t). Accordingly, the matrix elements of the 
transformed quasienergy operator 

{{a'm'\Q'\am)) = {a'\H'^,_^\a) d- Sm'mSa'amhuj, 


in reference 1511. 


(26) 
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Figure 2. By block diagonalization of the quasienergy operator with respect to 
the photon index m, one obtains the effective Hamiltonian Hp. 


are determined by the Fourier components H!^ = ^ dt H'(t) of the gauge- 

transformed Hamiltonian. 

The unitary operator Ud that diagonalizes the quasienergy operator with respect 
to a certain basis |am)), 

{{a'm'\Ul,QUD\arn)) = Sm'mSaiai{a\HD\a) mfuu) (27) 

is constructed such that it leads to a time-independent gauge-transformed Hamiltonian 

Hd = uUt)H{t)UDit) - thujj{t)dtUDit) (28) 

that is diagonal with respect to the basis states |a), 

{a'\HD\a)=d^,a.So.. (29) 

The Floquet Hamiltonian is related to Hd via the unitary transformation 

H[„=UDito)HDUj,ito). (30) 

The Floquet mode \uam)) = Hu jam)) with quasienergy £am = £« + rnhui reads 
I Warn (i)) = Hu(i)|c«)e™‘^‘, so that the micromotion operator can be expressed like 

UF{t,t') = UDmhit'). ( 31 ) 
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3. Block diagonalization of the quasienergy operator and effective 
Hamiltonian 


The quasienergy eigenvalue problem (201 is a convenient starting point for computing 


the Floquet Hamiltonian and the micromotion operator directly, without the need 
to compute the Floquet modes and their quasienergies. For this purpose one does 
not need to fully diagonalize the quasienergy operator. Instead one has to find a 
unitary operator Up that block diagonalizes the quasienergy operator with respect to 
the “photon” index to, 

{{a'm'\UpQUF\am)) = Sm'm{{a'\HF\a) + Sa'amHuj), (32) 

as illustrated in Fig.[^ Here we have introduced the gauge-transformed Hamiltonian 
Hf = uUt)H{t)UF{t) - ihUl{t)<ltUF{t), (33) 

which by construction is time independent. In fact, choosing an operator Up that block 
diagonalizes the quasienergy operator is equivalent to choosing Upit) such that the 
gauge transformation (33) leads to a time-independent Hamiltonian Hp- This time- 


independent Hamiltonian Hp is called effective Hamiltonian. Note that the unitary 
operator Up is not determined uniquely. For example, multiplying Uplt) with any 
time-independent unitary operator from the right leads to a mixing of states within 
the diagonal blocks of UpQUp, but does not destroy the block diagonal form. Unlike 
the operator Up that diagonalizes the quasienergy operator. Up does not depend on 
the basis states |a). 

Now, each of the diagonal blocks of UpQUp represents a possible choice for the 
Floquet Hamiltonian. This can be seen by writing the quasienergy operator like 
Q = Sm=-oo Wnm)){£a + mfuj) {{unm\ and Comparing it to the representation 
of the Floquet Hamiltonian. From the m = to' = 0 block one obtains 


H 


= ^ |a'0(<o))F i^((<a^0|(5|a0))F F(<a0(to)| 


(a|C/j)(to) 


(34) 


= Y.UF{to)\a') {{a'OplQUFH 

cx'a 

= UFUUjHpUpftf), 

where we have defined the rotated basis states 

\amj) F = UF\am'f) (35) 

and used Eq. (32). We can see that the Floquet Hamiltonian Hf^ is equivalent to 


the effective Hamiltonian Hp in the sense that both are related to each other by a 
unitary transformation. Moreover, we can use the unitary operator Up to construct 
the micromotion operator: 

UF{t,f) = Up{t)Ul{t'). (36) 

From Hf and UpU, t') one can then directly obtain the time evolution operator using 


Eq. (12). 


However, the time evolution operator U(t 2 ,ti) can also be express ed d irectly in 
terms of Hp and UpU) without introducing Ilf and UpU',t). Namely 42 , 

U(t2M) = (37) 

Compared to the representation (12) of the time-evolution operator in terms of the 
Eloquet Hamiltonian Hf and the micromotion operator f7F(^2,U), this expression 
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has the disadvantage that it is a product of three operators and not just of two. 


However, using the representation (371 has also advantages. The micromotion has been 


expressed by the one-point micromotion operator {7^(1), instead of by the two-point 
operator Upit, t'), and the phase evolution is described by an effective Hamiltonian Hp 
without the parametric dependence on the switching time to of Hf^^. The micromotion 
operator UF{t) can also be expressed like 

Upit) = exp (6(1)) (38) 

in terms of an anti-hermitian operator G = . The hermitian operator K{t) = iG(t) 


(39) 


has recently been given the intuitive name kick operator 43 
The diagonalization of the effective Hamiltonian Hp, 

provides the Floquet modes and their quasienergies: 

\Unmit)) = Up{t)\Un}e^^‘^\ (40) 

^nm — 4“ rnHu. (^1) 

Thus, the Floquet modes |u„(t)) = |m„o(^)); which describe the micromotion, are 
superpositions 

i»,.(*)) = E 'yan\o!(t))p. (42) 


of the time-dependent basis states 

|a(t))F = Up{t)\a0) = Up{t)\a), 
with time-independent coefficients 
"Tan (o:|rt^). 


(43) 

(44) 


The strategy of computing the effective Hamiltonian directly, without computing 
the Floquet states before, separates the Floquet problem into two distinct subproblems 
related to the short-time and the long-time dynamics, respectively. The first problem, 
computing the effective Hamiltonian (as well as the micromotion operator), concerns 
the short-time dynamics within one driving period only. The second problem consists 
in the integration of the time evolution generated by the effective Hamiltonian for a 
given initial state or even in the complete diagonalization of the effective Hamiltonian. 
This separation allows to address the long-time dynamics over several driving periods 
in a very efficient way, without the need to follow the details of the dynamics within 
every driving period. 

The advantage of splitting the Floquet problem into two parts becomes apparent 
especially when one of the two problems is more difficult than the other. A simple 
example for a case where computing the effective Hamiltonian is more difficult than 
diagonalizing it, is a periodically driven two-level system corresponding to a spin-1/2 
degree of freedom. While the block diagonalization of the quasienergy operator can 
generally not be accomplished analytically, the effective Hamiltonian describes (like 
every time-independent 2x2 Hamiltonian) a spin 1/2 in a constant magnetic field 
leading to a simple precession dynamics on the Bloch sphere. Thus, once the effective 
Hamiltonian and the micromotion operator are computed, the time evolution is known. 
An example for the opposite case, where the effective Hamiltonian can be computed 
at least approximately while its diagonalization is much harder, is a time-periodically 


driven Hubbard-type model 19 . It describes interacting particles on a tigh-binding 
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lattice. This driven model allows for a quantitative description of experiments with 
ultracold atoms in optical lattices. In the limit of high-frequency forcing a suitable 
analytical approximation to the effective Hamiltonian can be well justified on the time 
scale of a typical optical lattice experiment. However, the effective Hamiltonian will 
constitute a many-body problem that is difficult to solve. 

The possibility to compute the effective Hamiltonian for a many-body lattice 
system, at least within a suitable approximation, is also the basis for a novel and 
powerful type of quantum engineering. Here the properties of the effective Hamiltonian 
Hp are tailored by engineering the periodic time dependence of the Hamiltonian 
H[t). This Floquet engineering has recently been successfully applied to ultracold 
atomic quantum gases (see references in the introduction). The fact that the effective 
Hamiltonian can possess properties that are hard to achieve otherwise, like the 
coupling of the kinetics of charge-neutral atoms to a vector potential describing 
an (artificial) magnetic field 27 28 30 -32 34 -39 55 -60 , makes Floquet engineering 


also interesting for quantum simulation. Here, a quantum mechanical many-body 
model is realized accurately in the laboratory in order to investigate its properties by 
doing experiments. An essential prerequiste for Floquet engineering is an accurate 
approximation to the effective Hamiltonian. In the next section we will systematically 
derive a high-frequency approximation to both the effective Hamiltonian and the 
micromotion operator by block diagonalizing the quasienergy operator by means of 
degenerate perturbation theory. 


4. High-frequency expansion from degenerate perturbation theory 


Degenerate perturbation theory is a standard approximation scheme for the systematic 
block diagonalization of a hermitian operator into two subspaces—a subspace of special 
interest on the one hand and the rest of state space on the other—that are divided by 
a large spectral gap. Here we adapt the method such that it allows for a systematic 
block diagonalization of the quasienergy operator with respect to the “photon” index 
m (Appendix C). Moreover, we will identify the system-independent “photonic” 
—ihdt of the quasienergy operator (191, with {{a 


— ihdt\am}} 


part 


as the unperturbed problem. As a consequence the system-specific Hamiltonian H [t) 
constitutes the perturbation. This will allow us to systematically derive simple and 
universal expansions for both the effective Hamiltonian Hp and the micromotion 
operator Up{t) in the high-frequency limit, where hui constitutes a large spectral gap 
between the unperturbed subspaces (see Fig. [^. We would like to point out that the 
application of degenerate perturbation theory in the extended Floquet Hilbert space 
is a well established method. For example, it has recently been employed to estimate 
the matrix element for the resonant creation of collective excitations in a driven Bose- 
Hubbard model 51 and to treat a dissipative driven two-level system 61 . 

The basic strategy of our perturbative approach can be summarized as 
follows. The quasienergy operator is divided into an unperturbed part Qo and a 
perturbation V, 


Q — Qo + ^ ■ 


(45) 


The unperturbed operator can be diagonalized and separates the extended Floquet 
Hilbert space F into uncoupled subspaces Fm'^ of sharp “photon” numbers m with 
projectors Pm — |aw))((am|. These subspaces shall be separated by unperturbed 

spectral gaps of the order of Huj, which are assumed to be large compared to 





















High-frequency approximation for periodically driven quantum systems 


12 



Figure 3. Structure of the quasienergy operator Q = Qq + Vf) + = 

Qo + Id + the unperturbed problem is given by Qg^ ^ and the perturbation 
v(') can be separated into a block diagonal part ^ that conserves the “photon” 

number m and a part ^ comprising m 7 ^ 0-photon processes, (a) Generic choice 
of the unperturbed operator Qq. (b) Simple system-independent choice of the 
unperturbed problem Qo{t) = —ihdt to be used here; all unperturbed states of 
identical “photon” number m are degenerate. 


the strength p of the perturbation coupling states of different subspaces. When 
smoothly switching on the perturbation, such that the spectral gaps do not close, 
the unperturbed subspaces will be transformed adiabatically to the perturbed 
subspaces J-m corresponding to a diagonal block of the perturbed problem. Since 
the perturbation is weak compared to the gap, will differ from by small 
admixtures of states ^ only. This admixture will be calculated perturbatively 
by expanding a unitary operator Up that relates the basis states jam)) spanning 
the unperturbed subpsaces to the basis states \am))F = Uplam)) spanning 

the perturbed subspaces J-m- In contrast, if the spectral gap separating different 
subspaces would close, arbitrary weak coupling can hybridize degenerate states of 
different subspaces, contrary to the assumption of a weak perturbative admixture. 
The general formalism is developed in Appendix C| and will be applied to a specific 
choice of the unperturbed problem in the following. 

For the procedure described above a general and legitimate choice of the 
unperturbed problem would consist in the diagonal terms of the quasienergy operator 
with respect to a conveniently chosen set of basis states \a,m)), 


Qo = EE \am)) {{am\Q\am)) {{am\ 

m a. 

= EE |aTO))(e^°^ -b mHuj){{am\, 


m 


(46) 
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with Sa^ = {a\Ho\a). The operator Qq is diagonal with respect to the basis states 
jam)) by construction and the corresponding perturbation V' = Q — Qq consists of 
a block-diagonal part Vf, that couples states \am)) and \a'm)) of the same “photon” 
number m and a block-off-diagonal part that couples states \am)) and \a'm''fj of 
different “photon” numbers m' and m. The problem to be solved by perturbation 
theory is visualized in Figure [^a). The unperturbed problem and the perturbation 
expansion depend on the choice of the basis states |a). 

However, for the sake of simplicity we will not use Eq. (46). Instead we will 
simplify the unperturbed problem further, reducing it to the “photonic” part of the 
quasienergy operator, 


Qo{t) = -*Mt, 


(47) 


or 


Qo = EE \am))mtiujilam\, (48) 

m a 

which does not depend on the system’s Hamiltonian. For this choice the unperturbed 
quasienergies are degenerate within each subspace and read Sam = mhu. So Qo is 
diagonal not only with respect to a specihc set of basis states, but with respect to any 
set of basis states of the type \am)). The perturbation is given by the Hamiltonian, 

V{t) = H{t) (49) 


or 

V = H = \a'm')){a'\Hm'-m\a){{am\. (50) 

m'm ol' ct 

It can be decomposed like 

V = Vd + Vx. (51) 

Here the block-diagonal part Vd comprises the m' = m terms describing zero- “photon” 
processes determined by the time-averaged Hamiltonian, 

Vd = Ho, (52) 

Vd = EE \a'm'j}{a'\Ho\a)llam\. (53) 

m a'a. 

The block-off-diagonal part Vx describes Am-“photon” processes determined by the 
Fourier components Hxm of the Hamiltonian, 

Vx(t)= (54) 

Am^O 

= E E E |a'm -I- Am))(a'|i7Am|Q;)((Q!m|. (55) 

m Am^O cx'a 


The problem is visualized in Figure [^b). Its simple structure will allow us to write 
down universal analytical expressions for the leading terms of a perturbative high- 
frequency expansion of the effective Hamiltonian and the micromotion operator in 
powers of p/hw, with p symbolizing the perturbation strength. 

Before moving on, we note in passing that it can be useful to shift the “photon” 
number of an unperturbed state \am)) by some integer Am, before applying the 
high-frequency approximation. Such a procedure can be useful, if two states |ai) 
and \a2) have time-averaged energies Sa'^ = {a\Ho\a) that are separated by Am 
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“photon” energies oi Huj, so that — Sai = Amhuj + 6 with 6 hui. In this 
case the two unperturbed basis states \aim)) and |a 2 m)), which are degenerate 
with respect to the unperturbed quasienergy operator, have average quasienergies 
{{am\Q\amf} = + mhuj that are also separated by the large distance Amhio + 5. 

Obviously, this violates the requirement that the perturbation should be weak. In 
turn the states \aim)) and \a2{rn — Am))) have average quasienergies that are nearly 
degenerate. Thus it is useful to redefine the “photon” number of states with quantum 
number a 2 , so that \a2m)y = \a2{m—Am))). This redefinition is equivalent to a gauge 
transformation (25) in H, where the unitary operator e~'^^'^^*\a 2 ){a 2 \ is employed 
to shift the time-averaged energy of \a2) by —Amhw. After this transformation 
the high-frequency approximation can be applied and used to describe the resonant 
coupling between both states |ai) and |a 2 )- Such a procedure can be employed, e.g., 
to describe resonant “photon”-assisted (or AC-induced) tunneling against a strong 


( 0 ) 


potential gradient 12 45 


4-1. Micromotion 

We wish to compute the unitary operator Up that relates the unperturbed basis states 
jam)) to the perturbed basis states \am'fjF that block diagonalize the quasienergy 
operator in a perturbative fashion. In the canonical van Vleck degenerate perturbation 
theory, it is written like 

Up = exp(G'), (56) 

with anti-hermitian operator 

G = -GU (57) 

In order to minimize the mixing of unperturbed states belonging to the same 
unperturbed subspace, it is, moreover, required that G is block off diagonal. One 
can now systematically expand G like 




(58) 


in powers of the perturbation. The general formalism for the perturbative expansion of 
G in a situation where the state space is partitioned into more than just two subspaces 
is described in [Appendix C[ Differences with respect to the standard procedure, 
where the state space is just bipartitioned, arise as a consequence of the fact that 
for multipartitioning it is generally not true anymore that the product of two block- 
off-diagonal operators is block diagonal. 


The general form of the leading terms of the expansion (58) is given by Eqs. (C.40) 


and (C.41) of Appendix C Let us evaluate them for the particular choice of the 


unperturbed problem ( |47| ). Apart from 
((a'mlG^'^^lam)) = 0, 


(59) 

for all diagonal matrix elements, following directly from G being block-off-diagonal, 
for m' m we obtain 


((a'm'IG^^^Iam)) = — 


yCr \H jji'— 

im' — m)fiui 


(60) 
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and 


= (°'l l°> 

[(m' — m)hu\^ 


E 

1 


{(X \Hjyi> Hjji'' — jyi\(x) 

(m' — m)Hu! 

1 


(to" — m')huj (to" — m)hw 
We can now also expand the unitary operator 17^;’ in powers of the perturbation, 




{l') 

F • 




One finds 


uf = 1 , 


uP = g( 2) + i[G(i)]^ 


(61) 

(62) 

(63) 

(64) 

(65) 


( 66 ) 


where the second term of the last equation possesses matrix elements 
(a'm'|l[G<«]Vm))= 

2■' im' — — m){riUjY 

/ \ / 

which are finite also for m' = m. 

The corresponding operators in 1-L can be constructed by employing the relation 
A{t) = ^e*'"'“‘|a')((aW|A|aO))(a| (67) 

m 

that is valid for operators A that are translationally invariant with respect to the 
“photon” number, ((o'to + AtoIAIoto)) = ((aAml^laO)). In doing so. Up and G 
translate into time periodic operators Upit) and G(t) and Eq. (56) into 

Uplt) = exp(G(t)) (68) 


(see also Appendix B[ ) . The leading terms of the perturbation expansion take the form 

r.imujt 




m^O 


and 


G(2)(t) = ^ . 

m/0 


Cf (t) = i, 

i 7 y(t) = GW(t), 

r(2) 


e^mut [Ho, Hr, 


{mhwY 


\ E 


(69) 


2 m(m — m')(hw)‘^ 


uP{t) = G^^\t) + -Y, E 


„i{m+m')bjt TT TT 


2 2 ^ m'mihujY 

m^Om'^0 ' ' 


(70) 

(71) 

(72) 

(73) 
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One can express these terms also as time integrals. For the leading order we obtain 


G^^\t) = U^p\t) 


ntQ-\-T 


T 


to 


dt'H{t')Y, 


2i sm{muj{t — t')) 


m—1 


in 1 
hoj T 


rt+T 


dt 


'Hit')[ 


1 - 2 - 


mhuj 
J-t' 


T 


(74) 


In the final result we have separated a factor of representing the inverse integration 
time. It was obtained by setting the free parameter tg to tg = t allowing us to use 


E °° sin(fca;) tt — x ^ ^ 

for 0<x< 27r, (75) 

k=l 

which is formula 1.441-1 of reference [6^ . 

One can now approximate Up{f) up to a finite order v by simply truncating 
the perturbative expansion of Upit) like Upit) « J2'^=oUp\t). However, this 
approximation has the disadvantage that it does not preserve unitarity at any finite 
order i). In turn, truncating the expansion of G(t) leads to an approximation 


Upit) - exp (fi 

that gives rise to a unitary operator ll^p"^{t) for every finite v. 

The unitary two-point micromotion operator can be written like 

Upit^t') = UpmUt') = exp {F{t,t')) (77) 

with anti-hermitian operator F{t, t') = —F^{t, t'). Expanding F{t, t') in powers of the 
perturbation, 

OO 

F{t,t')=Y,F^’'Ht,t'), (78) 


and comparing the epxansion of exp (E(t,t')) in powers of the perturbation with that 
of exp (G(t)) exp ( — G{t ')), one can identify 

F(^)(t,t') = G^^\t) - G^'^\t'), (79) 

F^^\t,t') = G^^\t) - G^^\t') - ^[G(i)(t),G(i)(t')], (80) 

and so on. This gives the explicit expressions for the leading orders 


F^^\t,t') 


F‘'^\t,t') 



(81) 






m' ^0^771 


m{m — m')(hu}Y 


1 Amujt—im'ojt' \ tt T-T ~\ 

1 e 

2 mm'ihujY 


( 82 ) 
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An approximation preserving the unitarity of the micromotion operator reads 

UF(tX) ~ exp = Ufit^t'). (83) 

i/=i ^ 


Jj..2. Effective Hamiltonian 


In order to obtain the effective Floquet Hamiltonian from Eq. (34), we need to compute 
the matrix elements (32) for m = m' = 0, 

H^,^ = {a'\HF\a) = {{a'G\ulQUF\aO)) = Qo,a'a- (84) 

Expanding these matrix eleme n ts in p owers of the pertu rb ation, the lea ding terms 

Evaluating 


Appendix C 


Qg are given by Eqs. (C.50I, (C.51), (C.52|, and (C.53) of j 
these expressions for the unperturbed problem (47), we obtain the perturbative 
expansion for the effective Hamiltonian 


Hf = Y.H^ 


(^) 
F ’ 




with Hp'^ = Fja'a ^he leading terms are given by 


hP = 0 , 

hP = Ho, 


H 


(2) HmH —m 


= E 


mhoj 


(85) 


( 86 ) 

(87) 

( 88 ) 


H 


m^O 

(3) _ —m 5 [Ho, HP] 


= E 

mp^O 


2{mhuiY 


E 

m' pO^m 




(89) 

One can express these terms also in terms of time integrals. The leading order is given 
by the time-averaged Hamiltonian, 


H 


( 1 ) 


1 


F = ^ / <HH{t). 


(90) 


The first correction takes the form 

hp =^ 




/O 


/O 


1 

rp2 

27V 1 


dti 


E rnhcv 

^ TnpO 

rti —imuj{ti — t2) 

, [g«o.g(fa)] 

^ mpo 

rii 


^ dt,J^Ph(l-2*-^jp-)[Hih),H{h)], 


(91) 


where the sum over m has been evaluated using Eq. (75) and where we have separated 
a factor of l/(2r^) representing the inverse integration area. In f'th order the effective 
Hamiltonian is approximated by 


Hf^Y.^f^^hP. 


12^0 


(92) 
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The results obtained here via degenerate perturbation theory in the extended 
Floquet Hilbert space are equivalent to the high-frequency expansion derived in 
references 42 ■ 44 by different mean^ 


4-3. Role of the driving phase 


An important property of the approximation (92) to the effective Hamiltonian is that 

(93) 

(94) 


it is independent of the driving phase. Namely, a shift in time 
H{f) H\t) = H{t-t'), 

which leads to 


Hm H' — e 




Hr, 


does not alter the perturbation expansion of Hp, 


H 




^ = H^. 


(95) 

This is ensured by the structure of the perturbation theory, which restricts the 


products HmiHm 2 • • • that contribute to H^^ to those with mi-\-m 2 + - • • mi, = 0. 


As an immediate consequence, also the approximate quasienergy spectrum, obtained 
from the diagonalization ol Hp, does not acquire a spurious dependence on the driving 
phase. In this respect, the high-frequency approxima tion obta ined by truncating the 
high-frequency expansion of Hp at finite order, Eqs. (76) and (92), is consistent with 
Floquet theory. 

A time shift does, however, modify the terms of the unitary operator Up(t) in 
the expected way, 

ulf\t)^Uj^''\t) = ulf\t-t'), 


since 




(96) 


(97) 


4.4- Quasienergy spectrum and Floquet modes 

From the approximate Floquet Hamiltonian one can now compute the quasienergy 
spectrum and the Floquet modes by solving the eigenvalue problem 

i7M|h„)H (98) 

One obtains 

£n ~ (99) 

and 

\Unit)) « ^ "'W)- (100) 

r M 

Here we have allowed that the order p' of the approximate unitary operator Up (t) 
describing the micromotion can be different from the order v of the approximate 
Floquet Hamiltonian Hp\ which determines the Floquet spectrum and the dynamics 
on longer times. This corresponds to the approximation 

K « U^pHtoWp^ul^'^Hto) = ( 101 ) 

§ There is a slight discrepancy, however, concerning the third-order correction to the effective 
Hamiltonian. The second term of our expression ( |89[ l is different from the corresponding term in 
equation (C.IO) of reference [43], where a spurious factor of 2 is found. 
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to the Floquet Hamiltonian 

The reason why it is generally useful to choose v' independent of v is the following. 
In high-frequency approximation the time evolution from to to t is described by 

\m) « e—(102) 

n 

with Cn = (wn ’’^^(to)|'0(to))- The accuracy with which the expression 
Cn |wn ’‘'\t)) captures the true micromotion of the system does not depend on 
the time span (t — to) of the integration, simply because this expression is time 
periodic. In turn, with increasing integration time (t — to), the approximate phase 
factors 6“*®" (*-*o)/?i deviate more and more from their actual value 
Thus, the longer the time span t — to the better should be the approximation e„ « Sn ^, 
that is the larger should be v. In contrast, the order z/' can be chosen independently 
of (t - to). 


5. Relation to the Floquet-Magnus expansion 


In this section we relate the high-frequency expansion of Hp and G{t) to the Floquet- 
Magnus expansi on [4^ (see also 47p8|[^ ). A discussion of this issue can also be found 
in references 42 43 Recently, the Floquet-Magnus expansion has been employed 
frequently for the treatment of quantum Floquet systems. The starting point of the 


Floquet-Magnus expansion is the form (12) of the time evolution operator. 


U{t,to) = UF{t,to) exp ( - -{t-to)Hf^ 


= exp 


[Fit,to)) 


exp 




(103) 


Then both F{t,to) and are expanded in powers of the Fourier transform of 

the notation P{t) = [/^(t, 0), 


46 


47 


the Hamiltonian. Note that in references 
A(t) = F{t, 0), and F = —j^Hq is used, implicitly assuming to = 0. 

The Floquet-Magnus expansion of F{t,to) is reproduced by our expressions (81) 
and (82). The Floquet-Magnus expansion of can also be obtained within our 
formalism. Namely, expanding in powers of the perturbation 


gives 


^to — 2^ ^to 1 

(104) 

/V(l) 

(105) 

= X + + XX\to), 

= Hf + X{to)X + XXito) 

(106) 

+ XitXPXHtto), 



and so on. From these expressions one obtains 


_ TT 

— rio, 


H, 


to 

F{2) 

to 


= E 


m^O 


mhuj 






[Ho,H„ 


(107) 

(108) 
(109) 
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and, in a subsequent step, also 

f 


= - 


dtiH{ti), 


H, 


F{2) 
^0 


'to 
‘I'K 1 

= *‘i. *dWi).ffte)]. 




^0+^1 


( 110 ) 


( 111 ) 


those of references 46 47 


where we have again employed Eq. (751. For to = 0 these expressions correspond to 


Truncating the Floquet-Magnus expansion after the finite order D, the Floquet 
Hamiltonian is approximated like 






*0 


( 112 ) 


i/=i 


However, even though it is derived from a systematic expansion, this approximation 
is plagued by the following problem. For any finite order v > 2, the spectrum of the 
approximate Floquet Hamiltonian possesses an artifactual dependence on to, 

or equivalently on the driving phase. This is not consistent with the spectrum of the 
exact Floquet Hamiltonian Hf, which is independent of the driving phase. In second 


order, the to dependence enters with the second term of Eq. (109). Let us consider. 


for example, a periodic Hamiltonian with even time dependence, H{t) = H{—t), so 
that Hm = H-rn- In this case ' vanishes for to being an integer multiple of tt/w, 
while it is generally finite for other values of tg- Therefore, generally the Floquet 
Hamiltonians and obtained from the Floquet-Magnus approximation for 

times to yf tg are not related to each other by a unitary transformation, as it is the 
case for the exact Floquet Hamiltonian. 


The origin of this spurious tg dependence lies in the fact that the expansion (104) 


of the Floquet Hamiltonian implies also an expansion Upit) = 1-1- Up'^ Up\t) • • • of 
the unitary operator Upit). At any finite order, such an expansion does not preserve 


r(2) / 


FM[y 


unitarity and, thus, the spectrum of the approximate Floquet Hamiltonian 
deviates from the f/th-order spectrum obt ained by diagonalizing the approximate 


effective Hamiltonian Hp given by Eq. (92). 


This observation can be traced back further to the ansatz (103) for the time 


evolution operator. Bi-partitioning the time-evolution operator into two exponentials 


like in Eq. (103) does not allow for disentangling the phase evolution from the 


micromotion. This is different for the tri-partitioning ansatz 


Ult,to) = UF{t) exp ( --{t-to)H^ 


uWo) 


= exp 


(g( 0 ) 


exp 


-^{t-to)H^ 


exp (-G(<g)), (113) 


which underlies the pertur bativ e approach presented in the previous section. In 
the tri-partitioning ansatz (113), first Up{to) transforms the state into a “reference 


frame” where by construction no micromotion is present. Then the phase evolution 
is generated by the effective Hamiltonian, before at time t the state is finally rotated 
back to the original frame by UF{t). In contrast as it appears in the ansatz 
(103), carries also information about the micromotion. This fact is somewhat hidden. 


when the to dependence of the Floquet Hamiltonian is not written out explicitly like 
in Ref. 


47 , where tg = 0 is assumed. 
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However, since we know that the effective Hamiltonian Hp and the Floquet 
Hamiltonian possess the same spectrum, we also know that, when expanding 
both Hp and Hf^ in powers of the inverse frequency, also the spectra will coincide 
up to this order. This means that the tp-dependent second term of Eq. (109) will 
not cause changes of the spectrum within the second order (oc Instead this 

second term can contribute to the third-order correction of the quasienergy spectrum, 
together with the terms of . This argument generalizes to higher orders. 

Let us illustrate our reasoning using a simple example. A spin-1/2 system shall 
be described by the time-periodic Hamiltonian 


H{t) = aSx + bcos{ujt)Sy, 

with spin operators Si and Fourier components 

b 

2' 


Ho = aSx 


Hi = H_i = -Sy, 


(114) 


Hm = 0 for \m\ > 2. (115) 


According to equations (87) and (88), in second order the effective Hamiltonian is 
approximated by 

rH) 


Hp 


hP = H)f> + H'p , 


(116) 


with 


= aSx 


Hp^ = 0. 


(117) 


Here the second-order term Hf,' vanishes since [Hm, H-m] = 0. In second order the 
quasienergy spectrum is, thus, approximated by 


e± 


d2] 


= ±--ha. 


(118) 


In contrast, the second-order approximation of the Floquet Hamiltonian based on the 
Floquet-Magnus expansion. 




frFMl2] _ ^F(l) ^F(2) 


(119) 


does contain a second-order term. Namely, from Eqs. (108) and (109) one obtains 
tF{i) 


H, 


to 


— (^Sx 7 


- ^ sm{ojto)Sz. 


ab 


( 120 ) 


It leads to the approximation of the quasienergy spectrum 


e± « = ± ^ha\/1 -b [{b/tj) sin(a;to)] ^ 

= ±^ha(l-i-^[{b/uj)sm{ujto)p -( 121 ) 

We can now make several observations that illustrate the reasoning of the previous 
paragraphs. First, we can see that coincides with eP within the order of the 

approximation. Deviations that occur are proportional to while our second-order 
approximation should provide the correct terms up to the power uj~^. Second, despite 
the presence of a finite second-order term proportional to does 

not contain a correction cx uj~^. This is consistent with the fact that hP and, as a 
consequence, also elj_ do not contain a second-order term. Third, we can see that, 
unlike the exact quasienergy spectrum, the approximate spectrum depends on 

the time to and, thus, also on the driving phase. However, this dependence on to 
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(or the driving phase) occurs only in terms oc that are not reproduced correctly 
within the second-order approximation. In a third-order approximation, the spectrum 
will be captured correctly and be independent of the driving phase up to the power 


,-2 


As a further example, we will discuss the circularly driven hexagonal lattice in 
the next section. There, we will see that the spurious driving-phase dependence of 
the Floquet-Magnus expansion will, additionally, also induce a spurious breaking of 
the rotational symmetry of the quasienergy dispersion relation (Section 6.3). Thus, 
even a weak tg dependence can seemingly change the properties of the system in a 
fundamental way. Therefore, the Floquet-Magnus approximation should be used with 
care. The high-frequency approximation derived in the previous section does not 
suffer from this problem. 


6. Example: Circularly driven hexagonal lattice 

In this section, we will discuss an instructive example of the physics of particles 
hopping on a hexagonal lattice [see Fig. |^a)] subjected to a circular time-periodic 
force 

F{t) = —F[cos{ujt)ej; -I-sin(wt)ey]. (122) 

For a system of charged electrons such a force can be realized by applying circularly 
polarized light, whereas for a system of neutral particles (atoms in an optical lattice 
or photons in a wave guide) it can be achieved as an inertial force via circular lattice 


shaking 25 26 38 ,^. The driven hexagonal lattice is particularly interesting, as it 
is the prototype of a Floquet topological insulator 27,57 -60 . It was pointed out by 
Oka and Aoki 


27 


that a non-vanishing forcing strength F opens a topological gap 
in the band structure of the effective Hamiltonian. As a consequence, the system 
possesses a quantized Hall conductivity, when the lowest band is filled completely 
with fermions. While the original proposal [27| is considering graphene irradiated 
by circularly polarized light (see also Ref. [64|), the topologically non-trivial band 
structure described by the effective Hamiltonian has been probed experimentally in 
other systems: with classical light in a hexagonal lattice of wave guides 


ultracold fermionic atoms in a circularly shaken optical lattice 38 


39 and with 


We have decided to discuss the circularly driven hexagonal lattice here, even 
though its single-particle physics has been described in detail already elsewhere 
27 38 58 65 , because of several reasons. First, it is a paradigmatic example of a 


system where the second-order high-frequency correction to the effective Hamiltonian 
gives rise to qualitatively new physics. Second, since both directions, x and y, are 
driven with a phase lag of 7r/2, the model is suitable to illustrate the difference between 
the high-frequency expansion advertised here and the Floquet-Magnus expansion. And 
third, it allows us to set the stage for the ensuing discussion (see Sect. [7| of the role 
of interactions, which have been discussed controversially recently 66 6H. This issue 
includes two aspects: the impact of interactions on the validity of the high-frequency 
expansion as well as how interactions appear in the high-frequency expansion. 

Let us consider the driven tight-binding Hamiltonian 


Hdr{t) = - 


E 


J a^i 




Vi{t)hi. 


(123) 


The first term describes the tunneling kinetics, with the sum running over all directed 
links {£'£) connecting a site £ to its nearest neighbor £' on the hexagonal lattice depicted 


























High-frequency approximation for periodically driven quantum systems 


23 




Figure 4. Hexagonal lattice with sublattice A (blue) and B (red), (a) The 
Hamiltonian of the undriven model possesses real tunneling matrix elements 

— J between neighboring sites, (b) The effective Hamiltonian Hp Ri H^p'^ +H^p'^ of 
the driven system features modified real tunneling matrix elements —jH) between 
nearest neighbors, contained in H^p\ and complex tunneling matrix elements 

— (or —I |e“*®), contained in for tunneling in anticlockwise (or 

clockwise) direction around the hexagonal plaquette. 


in Fig. I^a). Here is the annihilation operator for a particle (boson or fermion) at 
the lattice site i located at and the tunneling parameter J is real and positive. 
The second sum runs over the lattice sites and describes the effect of the driving force 
in terms of the time-periodic on-site potential V(,{f) = —ri ■ F{t) and the number 
operator hg, = a\dg. The direction of the vector pointing from site £ to a neighbor £' 
defines an angle ipi'g, 

rg> -rt = a[cos[Lptie)e^ -I- sin((^£/^)ey], (124) 

with ipui = ipgig d-TT. This angle determines the temporal driving phase of the relative 
potential modulation between both sites, 

vii{t) — V£{t) = Facos{ojt — (fi'g). (125) 


6.1. Change of gauge 


As will be seen shortly, we are interested in the regime of strong forcing, where the 
amplitude K = Fa of the relative potential modulation between two neighboring sites 
is comparable to or larger than hoj. Therefore, the Hamiltonian H^rit) is not a suitable 
starting point for the high-frequency approximation. 

A remedy is provided by a gauge transformation with the time-periodic unitary 
operator 


25 


U (f) = exp 


(126) 

















24 


High-frequency approximation for periodically driven quantum systems 


where 




Frj 

fVjJ 


[— sin(a;t)ea; + cos(ujt)ey]. 


(127) 


This gauge transformation induces a time-dependent shift in quasimomentum, and 
the second integral has been included to eliminate an overall quasimomentum drift. It 
provides a constant that subtracts the zero-frequency component of the first integral, 
thus making the time average of xii^) over one driving period vanish. One arrives at 
the translationally invariant time-periodic Hamiltonian 


ff{t) = U'<{t)Hd,{t)Uit) - ihtj\t)U(t) = - ^ (128) 

{I'f.) 


Here the scalar potential vi{t) is absent while the driving force is captured by the 
time-periodic Peierls phases 

Ste{t) = xi{t) - XI' {t) = ^ sin(a;t - ipe'e). (129) 

Now we are in the position to apply the high-frequency approximation, even for 
K ^ hio. The actual requirement is that Huj must be large compared to the tunneling 
matrix element J, which determines both the spectral width of Hq and the strength 
of the coupling terms Hm with m Q. 


6.2. Effective Hamiltonian 


The leading term in the expansion of the effective Hamiltonian is according to Eq. (871 
given by the time-average of the driven Hamiltonian 

h''P =HQ = -Y,J^^^a\,d^. (130) 

{VI) 

It corresponds to the undriven Hamiltonian with a modified effective tunneling matrix 
element 

= ( 131 ) 

where Jn denotes a Bessel function of integer order n. This result was obtained by 
employing the relation 


exp(ir sin(s)) = 


E 

k— — oo 


Jk{r) exp(ifcs). 


(132) 


This Bessel-function-type renormalization of the tunnel matrix element, see Fig. [^for 
a plot, allows to effectively reduce or even “switch off” completely the nearest neighbor 
tunneling matrix element. This effect is known as dynamic localization [^, coherent 
destruction of tunneling [4][^ , or band collapse . It has been observed in the coherent 
expansion of a localized Bose condensate in a shaken optical lattice [^. The effect 
has also been used to induce the transition between a bosonic superfluid to a Mott 
insulator (and back) by shaking an optical lattice 19 The possibility to make 


the tunneling matrix element negative has moreover been exploited to achieve kinetic 
frustration in a circularly forced triangular lattice and to mimic antiferromagnetism 


with spinless bosons 25,26 
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The second-order contribution to the effective Hamiltonian is given by Eq. (88) 
and can be written like 


H 


( 2 ) 


oo ^ 

mhu} ^ 


m—1 




(133) 


(134) 


with the Fourier components of the Hamiltonian reading 

(t'l) 

By using the relation [d^d;, d)„d„] = 6im,a\a„ — Skndm^i ; which holds both for bosonic 
and fermionic operators d,, as well as J-ra{x) = one arrives at 


H 


( 2 ) 


d!,d« 


(135) 


where the sum runs over next-nearest neighbors and i. The effective tunneling 
matrix element is given by 


jZ 

= -'^m(^)2zsin(m[(^f.fc - ip£k]), 

m—1 


(136) 


where k denotes the intermediate lattice site between i' and f, via which the second- 

order tunneling process occur^ One can immediately see that the tunneling matrix 
( 2 ) 

elements purely imaginary and that they depend, as an odd function, on 

the relative angle pi'k — <ptk only. This relative angle is given by tpi'k — <Pik = 
with sign = -1-1 iotit = — 1) for tunneling in anticlockwise (clockwise) direction 
around a hexagonal lattice plaquette. Therefore, one finds 


(137) 

forming the pattern of effective tunneling matrix elements depicted in Fig. |^b). Here 
= -£ E (m2^/3) -^ 3 ^ ( 138 ) 


m—1 


and 


6 = sign(j(2^)|. 


(139) 


Since \ Jm{x)\ decays like with respect to the order to, for sufficiently small K/huj 
the sum is to good approximation exhausted by its first term, as is demonstrated also 
in Fig. 
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H, 


m 


H 


( 2 ) 


^ ’ ( 140 ) 

as it is depicted in Fig. Bb)> directly corresponds to the famous Haldane model 
(see also reference [^), being the prototype of a topological Chern insulator 70 71] . 
The next-nearest neighbor tunneling matrix elements open a gap between the two low- 
energy Bloch bands of the hexagonal lattice, such that the bands acquire topologically 
non-trivial properties of a Landau level characterized by non-zero integer Chern 
number 72 . As a consequence, the system features chiral edge states, as they 


II In other lattice geometries several two-step paths between I' and i can exist, in this case one has 
to sum over all of them. 
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Figure 5. Effective tunneling matrix elements (solid black line) and 

I(J'^Ihu)) (dashed blue line), as well as the leading term 
contributing to j( J'^ jhw) (dotted red line). 


have been observed experimentally with optical wave guides 39 


conductivity, as it has been observed with ultracold fermionic atoms 38 


and a finite Hall 
which is 


quantized for a completely filled lower band. The fact that circular forcing can induce 
such non-trivial properties to a hexagonal lattice has been pointed out in Ref. 27 . 
This is the first proposal for a Floquet-topological insulator (later proposals 
include Refs. 57 59 ). These systems can be defined as driven lattice systems with 
the effective Hamiltonian featuring new matrix elements that open topologically non¬ 
trivial gaps in the Floquet-Bloch band structure. Such new matrix elements appear in 
second (or higher) order of the high-frequency approximation that capture processes 
where a particle tunnels twice (or several times) during one driving period and that 
are of the order of ^ Therefore, Floquet topological insulators require the 

driving frequency to be at most moderately larger than the tunneling matrix element 
J. This is different for another class of schemes for the creation of artificial gauge fields 
and topological insulators recently pushed forward mainly in the context of ultracold 

In these schemes non-trivial effects enter already in the 


28-32 34-37 


quantum gases 

leading first-order term of the high-frequency expansion, such that they work also in 
the high-frequency limit hui ^ J. 

The non-interacting driven hexagonal lattice considered in this section can easily 
be solved numerically, without further approximation. The driven Hamiltonian (1281 
obeys discrete translational lattice, with two sublattice states per lattice cell. As a 
consequence, the single-particle state space is divided into uncoupled sectors of sharp 
quasimomentum, each containing two states. The problem is reduced to that of a 
family of driven two-level systems labeled by the quasimomentum wave vector k. The 
high-frequency approximation (140) is nevertheless useful. First, it gives rise to an 
analytical approximation to the Floquet Hamiltonian, directly corresponding to the 
paradigmatic Haldane model 38,58 . Second, as we will argue in the next paragraph. 


the second-order approximation captures already the essential physics of the non¬ 
interacting system in the regime of large frequencies. And third, it allows for taking 
into account also some of the effects related to interactions (see section 7.11. 

Let us briefly argue why the approximate effective Hamiltonian (140) captures 
the essential properties of the full one for the translationally invariant non-interacting 
system in the limit of large driving frequencies. The two-level Hamiltonian acting in 
the single-particle space of states with quasimomentum k is represented by a 2 x 2 
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matrix of the general form h{k, t) = /io(A;, t) + h{k^ t) ■ cr, with cr denoting the vector 
of Pauli matrices and vector h = /ly,/iz). If the elements hi{k,t) ^ J are small 
compared to huj, the perturbation expansion can be expected to converge. It is then 


left to argue that corrections beyond the second-order approximation (140) do not 
lead to qualitatively new behavior. 

Haldane 68 (see also reference 


This can be done employing the arguments by 
^). In the subspace of quasimomentum fc, the 


effective Hamiltonian Hp is represented by a time-independent matrix 
hp{k) = hFQ{k) -f hF{k) ■ cr. 

Its components determine the single-particle quasienergy dispersion relation 
e±{k) = hpoik) ± ^J\hF{k)\, 


(141) 


(142) 

The leading approximation of the effective 
H''p\ describes a hexagonal lattice with real nearest-neighbor tunneling 


with the two bands labeled by — and 
Hamiltonian 
matrix element 


J [see Fig. dlb)]. It gives rise to a contribution hp\k) to 


the vector hpik). Due to the fact that obeys space-inversion and time-reversal 
symmetry, two inequivalent quasimoment a fc + (the corners K and K' of the first 


Brillouin zone) exist, where h^^l{k±) = 0 69 . These are the Dirac points where the 


bands described by h^^\k) touch in a cone-like fashion. Moreover, also hi^\k) = 0 
for all k, since nearest-neighbor tunneling contributes only to the sublattice-mixing 
off-diagonal matrix elements. Now, the second-order correction to the effective 
Hamiltonian Hp contains complex next-nearest-neighbor tunneling matrix elements 
— [see Fig.&b)] of the order of J'^/{tvjj) [keeping K/{hbj) fixed]. This gives 


rise to a correction h}-p\k), where only the z-component h^pl{k) is non-zero, since 
next-nearest-neighbor tunneling is sublattice preserving. Moreover, the time-reversal 
symmetry breaking associated with the phase id makes this z-component non-zero 
at the Dirac points, such that h^p\{kp) = —h^p\{k_) 0. In this way the second- 

order correction removes the band touching and opens a gap between both bands. 
The opposite sign of hp^{kp) and hp^{k...) implies a finite Chern number of ±1 of 
the lowest band 


68 


Now it is important to note that in second order the bands 
are separated by a gap of order J'^/Ifiuj) (excluding situations, where the forcing 
strength is fine-tuned close to values giving = 0 or = 0). Even higher-order 
corrections, which in real space describe tunneling at even longer distances associated 
with a particle tunneling three or more times during one driving period, will be of the 
order of They will not be able to close this gap and change the topological 

properties of the bands anymore. Thus, the essential physics of the driven system is 


captured by the approximation (140). 


Qualitatively new behavior beyond the high-frequency approximation occurs, 
however, when fuo becomes small enough so that e+(fc) — e^{k) = mtuu with 
m = 1,2,3,... for some quasimomenta k. In this case both bands are coupled 
resonantly in an m-“photon” process and hybridize. Such a hybridization is not 
captured by the perturbative approach underlying the high-frequency approximation. 
It occurs, roughly, when the gap between subspaces of different “photon” number 
closes. The new band gaps resulting from the avoided quasienergy level crossing 
between the states of quasienergy s^{k) and e+(fe) — mHuj have been shown to give 


rise to intriguing topological properties without analog in non-driven systems 57,73 . 
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The circularly driven hexagonal lattice is also an instructive example that illustrates 
the difference between the high-frequency expansion of the effective Hamiltonian Hp 
on the one hand and of the Floquet Hamiltonian Hf^, as it appears in the Floquet- 
Magnus expansion, on the o ther. 


According to Eqs. (1081 and (1091, the leading terms of the expansion of Hf^ read 


H, 


HD ^ Ml) 


and 


H, 


F(2) _ M2) 
to — 


E 

m^O 


imujtQ 


mhu) 




(143) 


(144) 


Evaluating the difference between the Eloquet Hamiltonian and the effective 
Hamiltonian in second order, one obtains 


H, 


F(2) 

to 


- H 


( 2 ) _ 


= -E^. 

WD) 


( 2 ) ^ 
to{{t'l)) 


(145) 


Here 


J 


( 2 ) 




m^O 

huj 


( pirnMo-ipkl) _ 


ZJ 

- Ji(^) [cos(wto - Fki) - cos{ujto - (fitk)] 

(146) 

denotes the to"dependent part of the next-nearest-neighbor tunneling matrix element 

(147) 

It is easy to 


J 


( 2 ) 


+ 


( 2 ) 


{{£'£)) ^ ■^to{{Fl)) 

of with k defined as the intermediate site between £ and £'. 


check that the tunneling matrix element (1471 depends on the direction of tunneling 


and not only on whether a particle tunnels in clockwise or anticlockwise direction 
around a hexagonal plaquette. This directional dependence is determined by to and 
does not only concern the phase, but also the amplitude of next-nearest-neighbor 


tunneling matrix element (1471. The consequence is a spurious to"dependent breaking 


of the discrete rotational symmetry of the band structure of the approximate Floquet 


Hamiltonian -f This can be seen as follows. 

The driven Hamiltonian H {t) obeys the discrete translational symmetry of the 
hexagonal lattice so that quasimomentum is a conserved quantity. Therefore, and 
Hp do not only possess the same spectrum, but also the same single-particle dispersion 
relation e±{k). The symmetry of the hexagonal lattice with respect to discrete spatial 
rotations by 27r/3 is broken by the periodic force. However, the force leaves the 
Hamiltonian H{t) unaltered with respect to the joint operation of a rotation by 27r/3 
combined with a time shift by —r/3. This spatio-temporal symmetry ensures that 
the effective Hamiltonian Hp possesses again the full discrete rotational symmetry of 


F(2) 


the hexagonal lattice. This is reflected in the leading terms (130) and (1351 of the 


high-frequency expansion. As a consequence, the quasienergy band structure, given 
by the single particle dispersion relation e±(fc), is also symmetric with respect to a 
spatial rotation by 27r/3. The Floquet Hamiltonian whose parametric dependence 
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Figure 6. Quasienergy dispersion relation £+(fc) of the circularly driven 
hexagonal lattice, for hw/J = 4 and K/{hw) = 1.2. We compare the second- 
order high-frequency approximation of the effective Hamiltonian, [(a) and 

blue line in (c)] with the second-order Floquet-Magnus approximation, (fe), 

evaluated for uitg = 27r/3 [(b) and red line in (c)]. The Floquet-Magnus 
approximation (k) does not only acquire a spurious to dependence, but also 

breaks the discrete rotational symmetry of the exact dispersion relation £+(fc). In 
first order both approximations coincide, = £|j|^^(fc) [dashed line in (c)). 


on the time to indicates that it depends also on the micromotion, does not obey the 
discrete rotational symmetry of the hexagonal lattice when evaluated for a fixed time 
to- However, the single-particle spectrum of is still given by e±{k) and, thus, 
rotationally symmetric. This latter property of the exact Floquet Hamiltonian is 

not preserved by the Floquet-Magnus approximation Here 

the second-order term does not only lead to a spurious to dependence of the spectrum, 
but also breaks the rotational symmetry of the quasienergy band structure. This is 
illustrated in Fig. 6 where we compare the approximate dispersion relation e^^\k) 
resulting from the nigh-frequency approximation of the effective Hamiltonian with 
the spectrum obtained using the Floquet-Magnus approximation. In this 

respect, the approximate Floquet Hamiltonian is inconsistent with Floquet 

theory. The origin of this inconsistency is that the dispersion relation in Floquet- 
Magnus approximation depends on the driving phase, which, in turn, depends on the 
direction. Even though the spurious symmetry breaking should be small and of the 
order of J(J/ha;)^ [for {K/tujj) fixed], corresponding to the neglected third order, it 
still changes the property of the system in a fundamental way. Therefore, the Floquet- 
Magnus expansion has to be employed with care. 


6.4- Micromotion 


The micromotion operator Upit) = exp (G(t)) resulting from the periodic Hamiltonian 
H{t) is approximated by G{t) « -I- G^^\t). Here G^^\t) = 0 and, employing 


Eq. (69), we find 


G^^\t) 


E 


— p~ 


t TT 

-ti-r 


'^ge'e{t)dldg , 
{ti) 


mhu! 


(148) 
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with 


m—1 


^_ 'j'm 


2 J 

- - j-MD^os{ujt-^i,i). 


(149) 


Since ifu' = ipti + tt, it follows that gi'i{t) = —Puift), such that is anti- 

hermitian as required. With respect to the original frame of reference, where the 
system is described by the driven Hamiltonian i?dr(t), the micromotion operator is 
given by 


Ul{t) = UF{t)lI{t) « exp 



exp 




(150) 

The dynamics that is described by U{t) does not happen in real space, 
but corresponds to a global time-periodic oscillation in quasimomentum by 
^^[sin(wt)e 3 ; — cos(cLit)ey]. This momentum oscillation is significant when K hw 
and it is taken into account via the initial gauge transformation in a non-perturbative 
fashion, as can be seen from the Bessel-function-type dependence of the effective 
tunneling matrix elements on K/fnu. In turn, Upit) conserves quasimomentum 
and describes a micromotion in real space. This real-space micromotion becomes 
significant, when the tunneling time 2Tvh/ J is not too large compared to the driving 
period T, i.e. for J/hw not too small. A significant second-order correction Hf,' is a 
direct consequence of this real-space micromotion. In the high-frequency expansion, 
the real-space micromotion is taken into account perturbatively. 

By expanding the tunneling matrix elements also in powers of the driving strength 
K, one finds that 


J(2) ^ JK^ 

JW 


(151) 


The quadratic dependence on K indicates that large driving amplitudes K/huj ~ 1 
are important in order to achieve a topological band gap ~ | | that is signihcant 

with respect to the band width ~ At the same time, the linear dependence on 

the tunneling matrix element J reveals that moderate values of J/huj, for which the 
high-frequency expansion is still justified, can be sufficient (see also Fig. [^. 

The time-dependent basis states capturing the micromotion can be constructed 
from the basis of Fock states \{nt}) characterized by sharp on-site occupation numbers 
ng. They read 


\{ng}{t))F = UF{t)\{ne}) « exp 



\M) 


1 - 


E 

(e'e) 


2t7 \ 

(^) cos{ujt - ipFi)a\,dg j\{ni}), (152) 


For the second approximation we have expanded the exponential, such that the state 
is normalized only up to terms of order | J/hwp, and replaced gi'i{t) by the leading 


term of the sum (149). The micromotion is dominated by the oscillatory dynamics of 
particles from one lattice site to a neighboring one and back. The probability for a 
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particle participating in such an oscillation is of the order of | |^, with the 

coordination number z = 3 counting the nearest neighbors of a lattice site. 


7. Role of interactions 


If we consider a periodically driven quantum system of many particles, then the 
presence of interactions will influence the high-frequency expansion in two different 
ways. First, the interaction terms in the Hamiltonian will simply generate new terms in 
the high-frequency expansion. This effect will be discussed in the following subsection 
|7.1| for the example of the driven tight-binding lattice discussed in the previous section. 
Second, the presence of interactions will severely challenge the validity of the high- 
frequency expansion, since even if the single-particle spectrum is bounded with a width 
lower than hoj, this will not be the case anymore for collective excitations. This issue 
will be addressed in subsection El 


7.1. Interaction corrections within the high-frequency expansion 

Ignoring for the moment concerns that the high-frequency expansion should be 
employed with care in the presence of interactions, let us have a look at the new 
terms that will be generated by finite interaction terms in the Hamiltonian. We will 
focus on the example of the driven tight-binding lattice that was discussed on the 
single-particle level in the previous section. 

Starting from the interacting problem H^^if) -\- Hi^t with a time-independent 
interaction term Hi^t of the typical density-density type, the interactions are not 


altered by the gauge transformation (128). Thus, after the gauge transformation the 


interacting lattice system is described by the driven Hubbard Hamiltonian 

R(i) + Rint, 


with the time-periodic single-particle Hamiltonian H(t) given by Eq. (128). 
presence of the time-independent term Hint, with trivial Fourier components 


(153) 
The 

(154) 

will lead to additional terms in the high-frequency expansion of the effective 
Hamiltonian that we will denote by H^^^. The leading contribution appears in the 
first order [Eq. (87)] and is given by the time-independent operator Hint itself. 


Hi 


int m 






H 


( 1 ) 

Fint 


= H, 


The second-order correction (88) vanishes, 

H^fL = 0 , 


(155) 


(156) 


because all Fourier components Hntm with |m| > 0 vanish. Therefore, the leading 
correction involving the interactions appears in third order [Eq. (|89|)] and reads 


fV(3) _ 


OO 


771=1 


[H- 771; [Hint , .^m]] 

2{mhiL}Y 


h.c. 


(157) 


Here Hm denotes a Fourier component of the kinetic part H{t) of the Hamiltonian 
given by Eq. (134) of the previous section and “h.c.” stands for “hermitian conjugate”. 


The presence of interaction corrections beyond the leading order = Hint, as they 
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result from real-space micromotion, was overlooked in a recent work investigating the 
possibility of stabilizing a fractional-Chern-insulator-type many-body Floquet state 
with interacting fermions in the circularly driven hexagonal lattice 66 . A recent 


study, which is based on the results presented here, investigates the impact of the 


correction (1571 on the stability of both fermionic and bosonic Floquet fractional 
Chern insulators 74 . It is found that the correction tends to destabilize such a 


topologically ordered phase. 

Note that the high-frequency expansion of the Floquet Hamiltonian as it 
appears in the Floquet-Magnus expansion, produces also a second-order correction 


(109) that involves the interactions. It is given by 48 49 67 

„ 1 


H 


F{2) 


int,to 


= E 


imujto 


m—1 


mhuj 




(158) 


However, as we have discussed in sectionj^ if we approximate the Floquet Hamiltonian 
in second order this term will not influence the many-body spectrum within the order 
of the approximation, while it introduces an unphysical dependence of the spectrum 
on the driving phase. 

In order to get an idea of what type of terms appear in the high-frequency 
expansion, let us write down the third-order interaction correction for the simple case 
of spinless bosons in the circularly forced hexagonal lattice with on-site interactions 

Hint = ^y^,hiihe-l). (159) 


e 


This model system provides a quantitative description of ultracold bosonic atoms in 
an optical lattice and is interesting also because it might (as well as the fermionic 
system [^) be a possible candidate for a system stabilizing a Floquet fractional 
Chern insulator state. Using Eq. (157), we find those third-order correction terms 


that involve the interactions to be given by 
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(3) 

Fint 


= - ^ 2zWji^''hi{fie - 1 ) 


{te) 

+ E -ni- ht)d^ 


{t'M) 
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with coordination number z = i and coupling strengths 
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Figure 7. Like Fig. but with the spectral width defined by larger 

than the energetic separation hu) of sectors with different “photon” number m. 
The perturbation expansion is not justified and cannot be expected to converge. 
However, defining the unperturbed problem like in subfigure (b), the perturbation 
expansion can still be written down formally and, under certain conditions, still 
provide a useful approximation. 


(164) 

The third sum runs over all directed three-site strings (£'fc£), defined such that and 
k as well as k and I' are nearest neighbors, while t I'. The first sum produces 
a correction that reduces the on-site interactions of the leading term = H-mt- 

The second sum introduces both nearest-neighbor density-density interactions like in 
the extended Hubbard model and a pair-tunneling term. And the third sum, finally 
describes density-assisted tunneling between next-nearest neighbors, as well as the 
joint tunneling of two-particles into or away from a given site. 

7.2. On the validity of the high-frequency approximation for many-body systems 

The high-frequency expansion of the effective Hamiltonian Hp and the micromotion 
operator Uplt) can generally not be expected to converge for a system of many 
interacting particles. Namely, the time average Hq of the full many-body Hamiltonian 
H{t), which determines the spectrum of the diagonal blocks of the quasienergy 
operator Q depicted in Fig. [ij will possess collective excitations also at very large 
energies. Therefore the energy gaps of hu, which separate the subspaces of different 
“photon” numbers m in the unperturbed problem Qq, will close when the perturbation 
is switched on (unless hu was a macroscopic energy). 

Nevertheless, the fact that the unperturbed problem Qo is given by the “photonic” 
part of the quasienergy operator, with exactly degenerate eigenvalues mhu in each 
subspace of photon number m, allows one to formally write down the perturbation 
expansion. The energy denominators will not diverge. This is illustrated in Fig. 
And, even if the perturbation expansion can generally not be expected to converge, 
these terms can still provide an approximate description of the driven many-body 
system on a finite time scale. This can be the case if the i/th-order approximate 
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effective Hamiltonian is governed by energy scales that are small compared to huj. 
Then the creation of a collective excitation of energy hia corresponds to a significant 
change in the structure of the many-body wave function and is a process associated 
with a very small matrix element only. As a consequence, on time scales that are small 
compared to the inverse of such residual matrix elements, the approximate effective 
Hamiltonian can be employed to compute the dynamics and approximate Floquet 
states of the system. This is the basis for Floquet engineering in interacting systems. 

Let us illustrate this rather abstract reasoning by a concrete example. Consider 
again spinless bosons in the circularly forced hexagonal lattice, described by the driven 


Bose-Hubbard Hamiltonian (1531. The first-order approximation to the effective 
Hamiltonian is given by 

U 
~2 


(165) 


{te) 

In the thermodynamic limit (where the number of lattice sites M is taken to infinity, 
while keeping the filling of n particles per lattice site fixed) the system possesses 
excited states at arbitrarily large energies. Therefore, the width of the spectrum is 
obviously larger than Huj. As a consequence, the first-order quasienergy spectrum 
that describes the subspace of “photon” number m is much wider than Hu. It is not 
separated by a spectral gap from subspaces of different “photon number” anymore, but 
overlaps with quasienergies from these spectra. Therefore, the perturbation expansion, 
which is based on the assumption that the unperturbed spectral gaps do not close in 
response to the perturbation, can generally not be expected to converge]^ Namely, 
if we would consider the full diagonal blocks of the quasienergy matrix (Fig. as 
unperturbed problem, that is if we would add the block-diagonal part of the supposedly 
perturbation to the unperturbed quasienergy operator considered before, already 
the unperturbed quasienergy spectra of different subspaces m would overlap. As a 
consequence the energy denominators appearing in the perturbation expansion would 
diverge, whenever a degeneracy between unperturbed states of different “photon” 
number m occurs. 

We can now identify processes that spoil the high-frequency expansion for the 


driven Bose-Hubbard model (1531 and estimate the time scale on which they occur. 


For simplicity, we focus on the limit of strong interactions U ^ J and assume a 
filling of n = 1. In the spirit of the argument presented at the end of the preceding 
paragraph, we can now add part of the perturbation to the unperturbed problem. We 
include the interactions in the definition of the unperturbed quasienergy operator. 


Qq — Hint iHd-t^ 


so that the perturbation is given by the kinetic part (128) of the Hamiltonian, 

V'{t) = H{t)=Y,Hme^"""\ 


(166) 


(167) 


% Note that there are exceptions to this expectation. For example, the non-interacting system with 
U = 0 is described by a quadratic Hamiltonian. In this case the problem can be reduced to a 
single-particle problem. As a consequence, convergence is expected already for hu sufficiently large 
compared to the width of the single-particle spectrum ~ J. This is true even if hu is small compared 
to the width of the A^-particle spectrum ~ NJ of the Hamiltonian l |165| | with U = 0. Namely, the 
matrix elements for the joint creation of several single-particle excitations at once vanish. However, if 
the particles are coupled by a finite interaction strength U, the single-particle picture does not apply 
anymore. 
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with Fourier components Hm given by Eq. (134). The unperturbed Floquet states 
\{n(}m'fj are characterized by sharp site occupation numbers ni as well as sharp 
“photon” numbers m and read 

|{njm(t)) = e*"“^‘|{n£}). (168) 


They diagonalize the unperturbed problem, 

{{{n'i}m'\Q'Q\{nf}m)) = ( 169 ) 

where the unperturbed quasienergies read 

4nl}m = ^ H “ 1) + (1^6) 

I 


The high-frequency approximation is spoiled whenever two unperturbed states 
lln^jm)) and \{n'f}m''fj with different photon numbers m m' are nearly degenerate 
and coupled to each other (either directly or in a higher-order process). In such a 
situation a strong hybridization between these states of different “photon” number 
occurs, rather than a perturbative admixture of one state to the other as required by 
the high-frequency approximation. 

Let us investigate in how far the unperturbed “ground state” of the m = 0 
manifold, |V'o))) suffers from such detrimental coupling. This state corresponds to a 
Mott-insulator with one particle per site. 


IV'o)) = \{ni = 1}0)). 


(171) 


It is the only unperturbed to = 0 state without any doubly occupied site, so that 
the interaction energy vanishes. Its quasienergy reads = 0. It is the approximate 


ground state of the approximate effective Hamiltonians (165) in the regime U ^ J 75 


We can now systematically study the most relevant processes that couple \iPq)) to states 
of “photon” number to 0. For that purpose, we will follow the procedure applied 
to a one-dimensional chain in reference [51| . 

The state |'0o)) is coupled directly (in first-ordeiQ to unperturbed states 

|(^2^i)to)) = -^bl^b^^\{ni = 1}to)), (172) 


with one particle-hole excitation on neighboring sites £i and £ 2 , ™ photons, and 
unperturbed quasienergy = U-\-mhw. The resonance condition = £ 0 °^ 

is fulHlled for frequencies 


hio = - with TO = —1, —2, —3,.... (173) 

TO 

They are labeled by the negative change of the “photon” number, to, associated with 
the transition \ipo)) —>■ |(^ 2 ^i)w)). The corresponding coupling matrix elements is of 
the order of J and reads 




= 1 }) 

(174) 


The time scale for the resonant creation of a particle-hole excitations in the effective 
ground state can thus be estimated to be given by 2TrhjCm'^. If hiv is larger than U, 


■*" This order does not refer to the high-frequency approximation discussed so far, but to a perturbative 
approach treating |V^o}) collective excitations with m < 0 photons to be nearly degenerate. 
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the direct excitation of particle-hole pairs is not resonant, so that this type of heating 
process cannot spoil the high-frequency approximation. 

The next-strongest type of process is the creation of a collective excitation 
consiting of two coupled (i.e. overlapping) particle-hole pairs, with two extra particles 
on a site £3 and no particles on two neighboring sites £1 and £2 0 Such a state is given 

by 

|((£ 3 £i)(£'£ 2 ))m)) = ■^h\h\h^h^^\{ni = 1}to)) (175) 

and possesses the unperturbed quasienergy ~ -\-mtujj. The resonance 

condition ~ ^ 0 °^ fulfilled for frequencies 

with m = — 1 , — 2 , —4, —5, —7.... (176) 

Here we have omitted those m, for which already a dominating first-order resonance 
(173) occurs. The state (175) is coupled to \tpo)) only indirectly, in a second-order 
process via intermediate states |(£ 3 £i, 2 )w')) characterized a single particle-hole pair 
and photon number m!. The coupling matrix element can be evalutated using 
degenerate perturbation theory and reads 

E ((((^3£i)(£3£2))m|H(' E l(4£)m'))(((£3£)m'|)y|V-o)) 


m 


m — — 00 

1 

^ 2 

V6J2 


e=e,M 


1 


(0) _ (0) 

Jm-m' ilkij) 


+ ■ 


1 


(0) JO) 


e'o -e 


(^3^1,2)771' 


E 


hu ^ miZ — m’ 


(177) 


In order to arrive at a simple result, we have focused on the vicinity of the resonance 
and employed U = —mhu/Z in the energy denominators. The time scale for the 
resonant creation of such an overlapping pair of particle-hole excitations is given 
by 2nh/Cm . It is increased by a factor of the order of hio/J with respect to the 
creation of a single particle-hole pair, since it appears one order higher in perturbation 
theory. The experimentalists can avoid heating processes associated with the creation 
of two coupled particle-hole excitations by choosing huj well above 3£/. In this case, 
the most detrimental process will be the creation of a collective excitation given by 
three overlapping particle-hole pairs in a third-order process, with matrix element 

~ These third-order heating processes can, in turn, be avoided for huj 

well above 12U, and so on. Thus, by increasing the driving frequency, the time scale 
for detrimental heating due to the creation of collective excitations increases. 

The order of perturbation theory in which detrimental heating processes beyond 
the high-frequency approximation occur increases like a power law with the frequency 
huj, at least for the model and the parameters considered. This suggests that the time 


* The coupling matrix element for the creation of two or more independent (i.e. non-overlapping) 
particle-hole pairs, corresponding to an unconnected diagram, can be shown to vanish. Otherwise 
the overall rate for the creation of isolated particle-hole pairs would grow in a superlinear (i.e. 
superextensive) fashion with the system size. 
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scale for heating due to the creation of collective excitations increases exponentially 
with the driving frequency (the higher the energy of a collective excitation the more 
complex it is and the smaller will be the matrix element for its creation). Such a scaling 
is favorable for quantum engineering based on the high-frequency approximation. 
However, one has to note that the driving frequency cannot simply be increased to 
arbitrarily large values, because with increasing driving frequency another type of 
heating process will become more relevant. Namely, the resonant creation of high- 
energy single-particle states neglected in a low-energy model will become more and 
more relevant. In the driven lattice system these states belong to excited Bloch bands 
above a large energy gap A that are not taken into account in the tight binding 
description. If the resonance condition A « mhu) is fulfilled, these states can be 
populated in m-photon processes 76 -78 , the time scale of which tends to increase 


with m. Thus, Floquet engineering requires a window of suitable driving frequencies 
that are both large enough to suppress heating due to the creation of collective 
excitations and small enough to suppress heating due to interband transitions. In 
optical lattice systems, the existence of a window of suitable driving frequencies has 
been demonstrated in a number of experiments (see first paragraph of section [^. 

In the preceding paragraph the requirement “to suppress heating” should be read 
as “to suppress heating on the time scale of the experimental protocol”. Thus, the 
window of suitable frequencies depends also on the duration of the experiments, which 
in turn is determined by the effects to be studied. This dependence on the duration of 
the experiment is reflected also, in the response of a system to parameter variations. 
An important experimental protocol in the context of Floquet engineering is, for 
example, the smooth switching-on of the driving amplitude. Here the aim is to start 
from the ground state of the undriven model and to adiabatically prepare the ground 
state of the effective Hamiltonian obtained using the high-frequency approximation. In 
order to understand such parameter variations, one can employ the adiabatic principle 
for Floquet states 50 . In order to achieve the desired dynamics in a driven many- 


body system, it is not only required that the parameter variation is slow enough to 
ensure adiabatic following with respect to the high-frequency approximate effective 
Hamiltonian. At the same time, the parameter variation must also be sufficiently 
fast. Namely the matrix elements for the resonant creation of collective excitations 
discussed above will lead to tiny avoided crossings in the quasienergy spectrum between 
states of different “photon” number. These anticrossings are not captured by the 
high-frequency approximation and, in order to avoid heating, have to be passed 


diabatically 51 


The window of suitable driving frequencies is not necessarily determined by 
heating processes alone. Another limitation comes in, if the second-order term of the 
high-frequency approximation, Hf, , is relevant for the model system to be realized via 
Floquet engineering. The opening of a topological band gap in the circularly driven 
hexagonal lattice discussed in section [6^ is a prime example. This band gap appears 
in second order and depends inversely on the driving frequency uj~^ [if K/{fvjj) is held 
constant]. Thus if the band gap is required to be larger than some minimum value 
(determined, e.g., by interactions, temperature, or the rate of a parameter variation), 
this sets an upper limit for the driving frequency. The existence (and the identification) 
of system parameters for which a window of suitable driving frequencies exists is an 
important issue of Floquet engineering. 

The reasoning of this subsection suggests that, when dealing with a large system 
of interacting particles, on long time scales the high-frequency approximation can be 







High-frequency approximation for periodically driven quantum systems 


38 


expected to break down due to heating. In fact, the expected generic behavior of a 
driven many-body system in the thermodynamic limit is that it eventually approaches 
an infinite-temperature-like state in the long-time limit 
are conjectured to exist 


81-87 


79,80 . Though, exceptions 


including integrable systems and systems featuring 
many-body localization, where the size of the many- body state space is effectively 
reduced via the segmentation into different subspaces. In the sense of eigenstate 
thermalization 88 89 , the conjectured infinite-temperature behavior is attributed to 


the resonant (m-“photon”-type) coupling, and the resulting hybridization, of states 
with very different mean energies. However, irrespective of the long-time behavior, the 
high-frequency approximation might still provide an accurate description of a driven 
many-body system on a finite time scale. 


8. Summary 

We have used degenerate perturbation theory in the exended Floquet space to derive 
a high-frequency expansion of the effective Hamiltonian and the micromotion operator 
of periodically driven quantum systems. This approach provides an intuitive picture 
of the nature of the high-frequency approximation and its limitations. We have, 
moreover, related our approach to the Floquet-Magnus expansion. We have discussed 
that the latter is plagued not only by a spurious to dependence of the quasienergy 
spectrum, but that it can also violate further symmetries of the exact result, like 
the rotational symmetry of a Floquet band structure. Finally, we have addressed 
the possibly detrimental impact of interactions on the validity of the high-frequency 
approximation beyond a certain time scale. 
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Appendix A. Existence of Floquet states 


One can show that the eigenstates of the time-evolution operator over one period have 
the properties of Floquet states. The eigenstates |■^/'„(^)) of the time evolution operator 
from time t to time t -\-T fulfill 

U{t -\-T,t)\'lfri{t)) = an(t)\^l)n{t)). (A.l) 

Since U{t-\-T, t) is unitary, the eigenvalues o„(t) are phase factors, |a„(t)| = 1, and the 
eigenstates |'0n(t)) can be chosen to form a complete orthogonal basis. It is left to show 
that the time-evolved states \ifn{t’)) = tj{t',f)\'if{t)) are eigenstates of tj{T 1',t') 
with the same eigenvalue, i.e. that a„(t') = a„(t) = a„ = For that pu rpose 

one introduces 1 = U{t,t')U{t',t) on both sides of the eigenvalue equation (A.l) 
and multiplies the equation by U{t',t) from the left. Furthermore, the periodic time- 
dependence of the Hamiltonian has to be employed by using U{t',t) = U{t' -\-T,t-\-T). 
One arrives at 


U{t' T,t')\lljn{t')) = an{t)\fjn{t')), 


(A.2) 
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such that, indeed, a„(t') = an(t). Thus, one has 

U{t-\-T,t)\tpn{t)) = (A.3) 

with real quasienergy £„ not depending on the time t. The Floquet states \ipn{t)) 
fulfill 

\Mt + T)) = e-^^-^/^\fj„{t)) (A.4) 

and can be written like 

IV'n(O) = (A.5) 

where the time-periodic Floquet mode, 

\un{t)) = = |u„(t -h r)), (A.6) 

has been introduced. 


Appendix B. Relations between operators in H and T 


An important class of operators A acting in J- are those corresponding to operators 
A{t) = A{t -\- T) that act in H and possess a periodic and local time dependence. 
Here local in time means that the operator does neither involve a time derivative nor 
an integral over a finite time span. Let us briefly summarize some properties of such 
operators. 

(i) The operator A corresponding to Aif) possesses matrix elements 

{{a’m’\A\am)) = {a'\Am'-ni\a) (B.l) 


that are determined by the Fourier components 

Ara=lp 
^ do 


(B.2) 


of Aft). The fact that the matrix elements (B.l) depend on the difference m' — m 


show that time-periodic time-local (TPTL) operators Aft) correspond to operators A 
that are translation invariant with respect to the “photon” number m. Vice versa, for 
every operator A being translationally invariant with respect to m, we can construct 
a corresponding time periodic operator 


m 


(B.3) 


acting in % via 

Am = '^\a'){{a'm\A\atSfj{a\. (B.4) 

Oi' ,a. 


An example for such a TPTL operator is the Hamiltonian Hft), giving rise to matrix 
elements ((Q;'m'|.H|aTO)) = An example for an operator that is not 

time-local is the “photonic” part Qpft) = —ihdtl of the quasienergy operator. Here 
we have explicitly written out the unity operator 1 in "H, which we usually suppress. 
The corresponding operator Qp in T possesses matrix elements {{a'm'\Qp\am)) = 
Sm'mda'aiT^ltiv that are not translationally invariant with respect to the “photon” 
number. 
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(ii) A TPTL operator such as the Hamiltonian H{t) that for all times t is 
hermitian in H. corresponds to an operator H in J- that is translationally invariant 
with respect to the “photon” number m and hermitian, and vice versa, 

H\t) = H{t) ^ = (B.5) 

For example, one direction (“^”) can be shown as follows: {{a'm'\W\am'f) = 
{{am\H\a'm'))* = {a\Hrn-m'\ot')* = = {a'\Hm'-m\a) = 

Ifa'm'\H\arrif}. Here we have employed = —Hm'-m, which follows from 

Hit) = Em e^"^^*Hm^HHt). 

(iii) A multiplication of two TPTL operators in H directly corresponds to a 
multiplication in J-, 

Ait) = Bit)Cit) A = BC. (B.6) 

The proof is straightforward. This implies also that 

Ait) = f {Bit)) ^ A = f{B), (B.7) 

where the function / is defined via its Taylor expansion. 

(iv) When viewed as a constant function in time, the unity operator 1 in "H, 
with (a'llla) = 6a',a, directly corresponds to the unity operator 1 in B, with 
((a'TO'|I|am)) = Sm'm^a'a- That is 

Ait) = i A = l. (B.8) 

(v) A TPTL operator [/(t) that for all times t is unitary in H corresponds to an 
operator U in B that is translationally invariant with respect to the “photon” number 
m and unitary, and vice versa, 

U\t)Uit) = Uit)U'<it) = 1 um = UU'< = 1. (B.9) 

This is a direct consequence of (iii) and (iv). 


Appendix C. Degenerate perturbation theory in the extended Floquet 
Hilbert space 


Degenerate perturbation theory is an approximation scheme that allows for the 
systematic block diagonalization of a hermitian operator into two subspaces separated 
by a spectral gap. Here we apply the canonical van Vleck degenerate perturbation 
theory 90 to the quasienergy operator Q in the extended Floquet Hilbert space B and 


generalize it into an approximation scheme for the systematic block diagonalization of 
Q into multiple (more than two) subspaces separated by spectral gaps. The generalized 
formalism is found to contain additional terms that do not appear in the standard 
scheme for bipartioning. The procedure is closely related to the dressed-atom approach 
described in reference 53 and has previously also been developed for the concrete 


example of a driven two-level system 61 


Consider a Floquet system with quasienergy operator 


Q — Qo + 


A = l, 


(C.l) 


split into an uperturbed part Qo and a perturbation XV. The dimensionless parameter 
A shall eventually be set to one, and has been introduced to keep track of the order in 
which the perturbation appears. The eigenstates of the unperturbed quasienergy 
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operator the unperturbed Floquet modes, and their eigenvalues Sam, the 

unperturbed quasienergies, are known and fulfill 

Qo\am)) = e^°l\am)). (C.2) 

The index m separates the eigenstates into multiple subsets. The states within 
each subset m are labeled by the index a and span the unperturbed subspace 
related to m. The quasienergies of two subsets m and m' shall be separated by a 
quasienergy gap that is large compared to the matrix elements of the perturbation V. 
When the perturbation is switched on smoothly, without closing the spectral gaps, 
the unperturbed subspaces will be transformed adiabatically to the perturbed 
subspaces Tm, corresponding to a diagonal block of the perturbed problem. These 
subspaces Tm will be spanned by new basis states \am))B that deviate from the 
unperturbed states {lam))} by small perturbative admixtures of states from other 
unperturbed subspaces. The states \am))B are decoupled from states that are not 
in Tm, 

B{{a'm'\Q\am))B = 0 for m'm, (C.3) 

but generally they are no eigenstates of the perturbed quasienergy operator Q, so 
that B{{o:'m\Q\am))B can be finite also for a' a. The task to be accomplished 
by degenerate perturbation theory is to find systematic expansions for both the 
perturbed basis states \am))B, i.e. for the unitary operator U that relates them to 
the unperturbed basis states via 

\o:m))B = U\am)), (C.4) 

and the matrix elements b{{ o-'mlQlam ))b describing the physics within the 
subspace Tm- 

The projectors Pm into the unperturbed subspaces Tm'^ are defined by 

Pm = 'Y^\am)){{am\ (C.5) 

CX 

and obey 

= P (C.6) 

m 

with the unity operator 1. They can be used to decompose any operator A like 

A = Ad + Ax (C.7) 

into a block-diagonal part 

AD=Y.PmAPm (C.8) 

m 

and a block-off-diagonal part 

ix = H Pm'APm. (C.9) 

m 

^ m 

The product of two block-diagonal operators is again block-diagonal, 

AdBd = {AdBd)d i-e. {AdBd)x = 0, 


(C.IO) 
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and the product of a block-diagonal and a block-off-diagonal operator is block-off- 
diagonal, 

AdBx = {AdBx)x he. {AdBx)d=0, (C-H) 

AxBd = {AxBjj)x i.e. {AxBd)d=0- (C.12) 

If, like in the standard form of degenerate perturbation theory, the state space is 
bipartitioned only, one also finds AxBx = {AxBx)d, but this relation does not hold if 
the state space is partioned into more than two subspaces. Instead, for multipartioning 
one generally has 

AxBx = {AxBx)d + {AxBx)x- (C.13) 

The fact that the second term on the right hand side is finite will give rise to additional 
terms in the perturbation expansion that do not appear in the standard formalism. 

We wish to block diagonalize the full unperturbed quasienergy operator Q by 
means of a unitary operator U, such that 

WQU = Qo + W, (C.14) 

with block-diagonal operator 


W = Wd i.e. Wx = 0. 

It is convenient to split Eq. (C.14) into its block-diagonal part 

(C.15) 

[uHQo + Cd + Vx)U]d = Qo + W, 
and its block-off-diagonal part 

(C.16) 

[UHQo + Vd + Vx)U]x=0- 

(C.17) 

The unitary operator U defines the new basis states via Eq. (C.4|. 
in terms of an anti-hermitian operator G, 

It can be expressed 

C/ = exp(G), G=-G^ 

(C.18) 

However, U is not determined uniquely, unless one requires as the additional condition 

G = Gx i.e. Gd = 0 

(C.19) 

that keeps the mixing of states within each subpace small, 

the canonical van Vleck form of degenerate perturbation theory. 
We can expand W and U in powers n of the perturbation V, 

This ansatz defines 

OC 

W = ^ A"IT("\ 

n—0 

with 

(C.20) 

W(") = i.e. WP = 0 

and 

OO 

(C.21) 

C7 = 

(C.22) 


n=0 

The terms 17^"^ can be related to the terms in the perturbative expansion of G, 


g = Y, a”g(") 

n—1 


(C.23) 
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with 


Q(n) ^ Q(n) Q{n) ^ q 

(C.24) 

and 


Q{n) _ _ 

(C.25) 

One has 


= 1, 

(C.26) 

c/w = gw, 

(C.27) 


(C.28) 

[7(3) = i [G(i)g(2) + G(2)G(i)] + i [G(i)] ^ 

(C.29) 

etc.. Plugging these expressions into Egs. (|C.16|) and (|C.17|), we can 
determine W and G order by order. 

iteratively 

In zeroth order we find from Eg. (C.16) that 


ie(°) = 0, 

(C.30) 

while Eg. (C.17|) reduces to 0 = 0. For the next orders we obtain 


[G(i\Qo] 

(C.31) 

P^\Qo] = [yD,G«] +i[F^,GW]^ 

[G^^\Qo] = [Ez.,G(^)] +i[[Ex,GW],GW]^ 

(C.32) 

+ 5Fx.G<g,-i[[v;v,<;<'T,G''>]_, 

(C.33) 

etc. from Eg. (C.I6). And with that, Eg. (|C.I7) gives 


IE(i) = Vd, 

(C.34) 


(C.35) 

IE(3) = i [Ex, G(2)] ^ ^ [[Ex, G(1)], G(1)] 

(C.36) 


where, in order to obtain the expression for W^‘^\ we have used that 
[[yx,GW]x,G(i)]z3 = [[Vx,G(%G^^'>]d, since [[Vx,GW]jy,&^)]^ = 0. For the 
terms , Qo] the first deviation from the standard bipartioning perturbation theory 
appears in second order and is given by the second term on the right-hand side of 
(Eg. |C.3^ . Also the terms in the second line of Eg. (C.33) are new. In the expansion 
of W the first deviation occurs in the third order and is given by the last term 
of Eg. dele ). 

From the commutator [G^"\Qo] we can construct all matrix elements of G order 
by order. Since Gd = 0, we know that 


((a'TO|G^"^|am)) = 0. 


(C.37) 


The non-vanishing matrix elements of G*-"^ couple states jam)) and la'm')) with 
m' m and obey 

((a'm'IG^"^ lam)) = —((am|G^"^ |a'm'))*, 


(C.38) 
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following from being anti-hermitian. By using that 

((a'm'l [G(”\ Qo] I a™)) = ( 4 °i “ |G(”V™))> 

one can compute the matrix elements with m' m oi the leading terms: 


((a'm'IG^^^Iam)' 

((a'm'IG^^^Iam)' 


{{a'm'\Vx\oi'rn)) 


do) 

" ot'm' 


do) 


44 

(C.39) 

(C.40) 


= E 


E 


'(a'm' I Vd I a"m')) {{a"m' \ Vx \ am)) 
fe(o) _ .(0) 1 G(o) _ e(o) f 

a'm' I Vx \a"m)) {{a"m\ Vd \ am)] 


E E 


Va'm' ^a"m) 

a'm' I Vx I a"m")) {{a"m" \ Vx \am)) 


JO) _ JO) 


1 

^ 2 


^ m, ra 

1 


do) 


— e 


(0) 


do) 


^a.m 


(C.41) 


In Eq. (C.41) one can explicitly see that the last term would not appear for a 
bipartition of the state space, since m" must be different both from m and m'. 

When approximating the unitary operator t/ up to a finite order n we have two 
possibilities, either 

(C.42) 


or 


U ~ exp (g'i) + G(2) + ... + g(")) = i/I” 

c7 ~ 1 + U^^'> + + • • • + 


(C.43) 

Both approximations coincide up to order n in the perturbation. However, the first 
approximation has the advantage that it preserves the unitarity, C/I”! is a unitary 
matrix also for finite n. In turn, the second approximation does not preserve unitarity 
for finite order. However, sometimes unitarity is not relevant and it is convenient 
to employ the second approximation to compute corrections to the perturbed basis 
states (C.4) of the form 


I amj) B = A" I amj) 

n=0 


(n) 
B ’ 


\am))^^'^ = 


am) 


Therefore, let us evaluate also the matrix elements of the leading terms : 
{{a'm'\U^°^\am)) = {{a'm'\am)) = Sa'a^m'm 


(C.44) 

(C.45) 

(C.46) 




E 


'(a'TO'|yD|Q;"m')^ 


'(a"TO'|Vx|Q!m)' 


G(0) _ .(0) ') (JO) _ JO) \ 

Va'm' t-anijt.am) 

^ {{a'm'\Vx\a"m)){{a"m\VD\am)) 

^ (JO) _J0l)fJ0) _jo) ) 

Oi" y^OL'm' ^am) \^cx^rn' ^oi"m) 

a'm' I Vx I a"m")) {{a"m" \ Vx \ am)] 


E E 


(JO) -JO)) (JO) _ (0) X 

V-a'm' a" m") 


M\tJo) 


(0) 


^ m, m' 


(C.47) 
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It is left to determine the matrix elements of the diagonal blocks of the quasienergy 
operator that describe the physics within the subspace Hm, namely 

Qm,a'a = B {{o:'m\Q\am)) B = {{a'm\U^QU\am)) 

= {{a'm\{Qo + I^)|Q;m)). 

Expanding them in powers of the perturbation, 




(C.48) 

(C.49) 


n=0 


we find the leading orders to be given by 

<9m!a'a = {{a'm\Qo\a)) = (5a',a4m, 

<9m!a'a = {{a'm\VD\am)), 

Qm]a'a = {{a'm\Vx\a"m')) {{a"m"\Vx\am)) 


(C.50) 

(C.51) 


cx" m' 
^ m 


1 

^ 2 


1 


( 0 ) _ ( 0 ) 

^ ol' m ^ a" m' 


^(0) ^(0) 
tarn 


(C.52) 


^ E E 

ol"ol"' rn' \ 

7^ m 

i{a'm\Vx\ a"m')) {{a"m' \Vd\ a"'m')) \ Vx \ am)] 


(0) ^ ,'^(0) _ ^(0) N, 

f ^ocm f 


({a'm\Vx\ a"m')) {{a"m' \ Vd \ a"'m')) {{a"'m' \ Vx \ am)) 

\ a"m' cx'm/ \ oc'"m' oc'm/ 

((a'm|i^ |a"TO')) {{a"m'\Vx\a"'m)){{a"'m\VD\am)) 
Va"m' Va"m' t.am) 

{{a'm\VD\a"m)){{a"m\Vx\a"'m')){{a"'m'\Vx\am)) \ 

fe(o) dfe(o) I 

\ am a m/ \ am am) J 

+ i^E E E I 

a"a"' rn m" \ 

7 ^ m ^ m, m' 

la'm\Vx I a"m')) {{a"m' \ Vx | a'"m")) ((a'"m'' | Vx\am)) 

3 




J0)^)(^(0) 


do) 


/ (0) _ (U) N ( W) 

Va"m' 


( 0 ) \ / ( 0 ) 


— £ 


(0) 


IpW -£^"^ 

\^a'm ^a'"m'') \^a"m' ^a"'m" 

3 

(M _M )(jo) _Jo) N 

\^a'm ^a"'m") \^a'm ^a"m') 
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